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Abstract An overview of the plasma environment near the earth is provided. We describe 
how the near-earth plasma is formed, including photo-ionization from solar photons and 
impact ionization at high latitudes from energetic particles. We review the fundamental 
characteristics of the earth’s plasma environment, with emphasis on the ionosphere and its 
interactions with the extended neutral atmosphere. Important processes that control iono- 
spheric physics at low, middle, and high latitudes are discussed. The general dynamics and 
morphology of the ionized gas at mid- and low-latitudes are described including electro- 
dynamic contributions from wind-driven dynamos, tides, and planetary-scale waves. The 
unique properties of the near-earth plasma and its associated currents at high latitudes are 
shown to depend on precipitating auroral charged particles and strong electric fields which 
map earthward from the magnetosphere. The upper atmosphere is shown to have profound 
effects on the transfer of energy and momentum between the high-latitude plasma and the 
neutral constituents. The article concludes with a discussion of how the near-earth plasma 
responds to magnetic storms associated with solar disturbances. 
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1 Introduction 

In this article we discuss the near-earth plasma environment, for which the ionosphere is 
the natural focus. Because this ionized gas co-exists with the neutral atmospheric gas to 
which it is coupled, the ionized/neutral gas system is more properly considered a partially 
ionized gas rather than a medium consisting of two distinct, independent fluids. We discuss 
the origins of the ionized component and describe its general morphology and dynamics, 
both in the framework of the earth’s magnetic field as well as with respect to its coupling 
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Magnetosphere 



Fig. 1 The various paths of solar energy to the earth’s upper atmosphere and ionosphere as well as its 
coupling to the magnetosphere above and the mesosphere and troposphere below 


to the neutral atmosphere and its motions. Because of the emphasis of this monograph on 
atmosphere-plasma coupling, and in order to keep the article manageable, we do not include 
a detailed discussion of the earth’s magnetosphere plasma, except in its role of exchanging 
energy and momentum with the partially ionized atmosphere at its base. 

We discuss two main mechanisms by which ambient or “thermal” plasma is created in 
the near-earth environment, each of which depends on an interaction of an external energy 
source with the earth’s extended atmosphere: (1) photo-ionization of the neutral gas via 
impinging extreme ultraviolet (EUV) solar radiation and (2) impact ionization of the neutral 
gas via energetic particle precipitation. A second form of impact ionization, due to meteor 
ablation, is also briefly discussed. Whereas photo-ionization is responsible for the creation 
of the primary part of earth’s ionosphere, particularly at mid- and low-latitudes, at high 
latitudes impact ionization by precipitating, energetic particles creates thermal plasma with 
different properties and features. Although their ionizing mechanisms are driven by different 
energy sources, the resulting plasma in each case is a consequence of energy or momentum 
impact with the upper atmosphere itself. Indeed, as we shall see, it is precisely because the 
earth has an extended atmosphere that it also has an ionosphere. In fact, any planet with a 
sensible, gravitationally-bound atmosphere and an ionizing source of radiation of sufficient 
energy and intensity, will consequently form an ionosphere. For planets which also have 
an intrinsic magnetic field and a substantial rotation rate, such as the earth, the resulting 
plasma environment may be expected to be particularly dynamic and rich in its complexity 
and physical processes. 

In order to illustrate the context of the near-earth plasma within the larger sun-earth 
system, we depict the main pathways by which energy and momentum are deposited and 
exchanged within the earth’s upper atmosphere in Fig. 1. Here, one path shows how the sun 
directly influences the upper atmosphere via EUV and X-ray radiation, both via heating of 
the neutral atmosphere, which results in its expansion, and also via photo-ionization of a 
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small fraction of this same neutral gas which thereby re-generates the ionosphere through- 
out the day. A second path shows how solar wind processes, including plasma flow with 
embedded interplanetary magnetic fields, interact with the magnetosphere to drive large- 
scale electric fields, currents, and energetic particles that subsequently interact with the 
underlying upper atmosphere/ionosphere. Such solar wind/magnetospheric input impacts 
the near-earth plasma primarily at high latitudes, although it may dramatically influence 
all latitudes during periods of magnetic storms. Although the globally averaged energy in- 
put by the solar wind/magnetosphere is not as high as the total solar photon irradiance, its 
impact in the upper atmosphere is nonetheless quite significant, and locally it can be the 
dominant energy source. The diagram in Fig. 1 also includes a sketch showing that the up- 
per atmosphere is also influenced by energy and momentum sources from below, namely 
from electrical discharges in the troposphere (lightning) and mesosphere (e.g., sprites), as 
well as via upward-propagating tides, planetary waves, and gravity waves. Such sources 
have profound effects on the characteristics of the mid- and low-latitude ionosphere in 
particular. The energy source that heats and sets in motion the lower and middle atmo- 
sphere is, again, solar radiation, but in this case, its ultraviolet, visible, and infrared wave- 
lengths. 

An outline of this article is as follows: We begin with a description of how the earth’s 
ionosphere is formed via photo-ionization, and its daily variation. We then discuss how 
these ionized particles are influenced by the earth’s magnetic field and how they couple to 
the neutral atmospheric gas via collisions, as represented by the mobility tensor for electrons 
and ions. The dynamics and morphology of the ionosphere at low and mid-latitudes are then 
discussed, followed by a separate section on the earth’s high latitude ionosphere including 
the aurora and the exchange of energy and momentum with the magnetosphere. Finally, we 
discuss the effects of magnetic storms on the ionosphere and conclude with a brief summary. 

There are several excellent references that describe the earth’s ionosphere and upper 
atmosphere. These include: Hargreaves (1992), Kelley (2009), Prolss (2004), Rishbeth and 
Garriott (1969), and Schunk and Nagy (2009), among others. 


2 Creation of the Earth’s Ionosphere via Photo-ionization and Its Daily Variation 

In this section, we describe how the ionosphere is created by photo-ionization associated 
with solar EUV and X-ray radiation that impacts the earth’s tenuous upper atmosphere. We 
also describe how some positive ions thus produced may transfer their charge to different 
neutral species and how others subsequently recombine with free electrons to form neutral 
gas particles. We conclude with a general description of how the mid- and low-latitude 
ionospheric density and temperature vary with local time and altitude. 

2.1 The Earth’s Extended Neutral Atmosphere 

In the Introduction, we discussed how the external sources of energy and momentum im- 
pinge on the earth’s atmosphere to create thermal plasma via two main mechanisms: solar 
radiation induced photo-ionization and impact ionization from energetic particles and mete- 
ors. Before discussing these ionization processes, we first review some of the fundamental 
features of the Earth’s upper neutral atmosphere in order to illustrate how they influence the 
resulting plasma parameters formed by these ionizing mechanisms. 

Whereas the earth’s atmosphere near the surface consists primarily of nitrogen molecules 
(N 2 ) and oxygen molecules (CF), the atmospheric constituents change markedly at higher 
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Fig. 2 Model distributions of neutral and ionic constituents as a function of altitude for a mid-latitude loca- 
tion (Wallops Island, Va., USA) at 14 hr Solar Local Time. Ion distributions below 80 km are not provided 
by the model 


altitudes. At about 100 km and above, the atmosphere ceases to be homogeneous and con- 
verts to free-molecular flow with a large mean-free path of the particles. Additionally, the 
atmospheric constituents above about 100 km directly absorb solar EUV radiation. Because 
of this absorption and the fact that these gases are very ineffective infrared radiators, the 
neutral temperatures increase significantly with altitude above 100 km. Consequently, this 
region of the upper atmosphere is often referred to as the thermosphere. 

The main process by which the upper atmosphere absorbs EUV energy from the sun is 
photo-dissociation. Diatomic oxygen is the primary molecule in the upper atmosphere that 
dissociates. This occurs when it is struck by photons with wavelengths of 242 nm or less, as 
expressed by the reaction: 

Ot + hv(k < 242 nm) -^0 + 0 + energy 

The creation of atomic oxygen is important because it is readily photo-ionized, as discussed 
later. 

The various species of atoms and molecules that make up the atmosphere are shown 
in Fig. 2, for mid-latitude, daytime conditions with moderate solar activity, as computed 
from the NRLMSISE-2000 model (Picone et al. 2002). The concentrations per unit volume 
of the atmospheric constituents and their variation with altitude affect the ionosphere in 
two ways: First, they determine the amount of neutral gas available for photo-ionization 
at different altitudes, and second, they determine to what depth the solar EUV and X-ray 
radiation might penetrate, as discussed in the next section. For example, a consequence of a 
decreased neutral atmospheric density is that photons can penetrate to lower altitudes. 

Both the neutral and plasma densities in the earth’s upper atmosphere vary considerably 
because of EUV solar radiation variations associated with solar activity. Unlike the main 
portion of the sun’s emission spectrum, for example where the visible light is emitted, the 
EUV portion of the spectrum originates at the base of the highly dynamic corona, where 
it undergoes large variations in accordance with solar storms and solar activity. In fact, the 
EUV radiation typically changes by factors of 2-3 (sometimes larger) within the solar cycle. 
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Fig. 3 Representative plasma 
density profiles for day and night 
at solar minimum and maximum 
conditions at a mid-latitude 
location (Hargreaves 1992). The 
colored arrows show the 
approximate penetration depths 
of the X-ray, EUV, and UV 
radiation from the sun 
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The EUV radiation is doubly important for determining the characteristics of the ionosphere 
as it: (1) is absorbed by the upper atmosphere and thereby heats the atmosphere, causing it 
expand, as discussed above, and (2) photo-ionizes a small fraction of the gas, creating the 
ionosphere, as discussed below. 

2.2 Photo-Ionization, Charge Exchange, and Recombination 

The main mechanism that creates the ionosphere is the photo-ionization of a tiny fraction of 
the earth’s extended, tenuous neutral upper atmosphere during the daytime. During the night, 
there is a very small amount of radiation input resulting from the high altitude scattering of 
the solar radiation into the shadow of the planet and other minor sources. Incoming solar 
EUV and X-ray radiation penetrate the upper atmosphere with sufficient energy to detach 
an electron from a neutral atom to create an ion and electron pair, as depicted here for 
oxygen, the main atmospheric gas ionized in this manner: 

Of/iv — > 0 + + e“ 

This radiation penetrates to different altitudes based on wavelength. Enhanced plasma layers 
are formed where the intensity of the incoming radiation, which falls off with decreasing 
altitude, and the amount of a particular gas constituent available to be ionized, which falls 
off with increasing altitude, are optimal, as discussed by Chapman (1931). For the earth, two 
main layers of ionization are created in the daytime ionosphere due to the photo-ionization 
of oxygen atoms: one formed by soft X-rays around 110 km, denoted as the E region and 
one formed by EUV radiation near 150 km denoted as the FI region. 

The resulting ionization plotted as a function of altitude is shown in the left hand side 
of Fig. 3, for a mid-latitude location (Hargreaves 1992). Note the strong dependence of the 
resulting plasma density on solar activity, which directly affects the strength of available 
EUV and X-ray radiation as discussed above, as well as how the profiles differ between day 
and night, particularly at the lower altitudes (discussed below). The incoming ionizing radi- 
ation is represented schematically on the right hand side of the figure. The actual penetration 
altitudes at a given latitude vary during the day with the solar zenith angle. For large solar 
zenith angles, grazing incidence corrections must be taken into account. (See discussion in 
Rishbeth and Garriott 1969 and references therein.) 
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The number density distributions of the different ion species which constitute the iono- 
sphere are shown in Fig. 2, taken from the International Reference Ionosphere (IRI) empiri- 
cal model (Bilitza and Reinisch 2008; Bilitza et al. 2011), for the same mid-latitude, daytime 
conditions used for the neutral density profiles shown in this figure. The curves show how 
the various ion constituents change with altitude. Indeed, the distribution of ions within the 
ionosphere is controlled to a large extent by the distributions of the neutral gas constituents. 
When comparing the ionized density profile with that of the neutral density profile shown in 
Fig. 2, notice that, for altitudes below about 300 km, the plasma number densities are 10 4 to 
10 6 times smaller than those of the neutral parameters. 

It is important to emphasize that photo-ionization is sometimes only the first step in a 
chain of processes by which the ambient ion gases are ultimately formed. For example, 
oxygen atoms are readily photo-ionized while atmospheric molecules, such as Oi and NO, 
do not photo-ionize nearly as efficiently. However, charge exchange readily occurs between 
the oxygen ions and the ambient nitrogen and oxygen molecules (No and Oi), creating NO + 
and Ot : 

0 + + N 2 NO + + N 
0+ + Oo -> ot + o 

At the lower altitudes, the resulting molecular ions are the dominant ones, even though the 
initial ionization is produced by photo-ionization of oxygen atoms. This is an important 
aspect of the near-earth plasma because the molecular ionization is lost quickly through 
recombination when their ionization source is cut off. 

As soon as the positive ions and electrons are freed, they are subject to chemical recom- 
bination processes as well as charge exchange and transport processes. The most common of 
these processes is dissociative recombination, which is important for the E-region molecular 
ions: 


Ot T e — > O T O 
NO + + e -* N + O 

This is a rather fast process in which newly formed 0| and NO + typically recombine in 
seconds. In fact, the loss rate for the molecular species, L, is given by: L = C r NjN e , where 
C r is the recombination coefficient, N, is the ion density, and N e is the electron density. 
Since the loss rate depends on the ion density, if the molecular species is the dominant ion, 
the loss rate is very rapid. On the other hand, 0 + recombines via a radiative process: 

0 + + e -> O + hv 

which is very slow. In fact, 0 + is far more likely to be lost via charge exchange or other 
processes rather than through recombination. Since the molecular ions dominate the lower 
altitudes (i.e., less than 200 km), when the direct illumination of the sun vanishes at dusk, the 
lower ionosphere, which rotates with the planet, becomes largely depleted of any significant 
plasma until local sunrise at which time the main photo-ionization process starts over again. 
In this fashion, we see that at the lower altitudes, the ionosphere is continuously being re- 
born every day. 

The physical processes and chemistry associated with the ion and neutral gas con- 
stituents, with particular emphasis on ionization and recombination, constitutes the general 
domain of aeronomy, the science of the upper atmosphere. A detailed discussion of these 
processes may be found in Rishbeth and Garriott (1969), Banks and Kockarts (1973), Prolss 
(2004), and numerous other texts. 
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PLASMASPHERE 



Fig. 4 The plasmasphere, an extension of the ionosphere at mid and low latitudes. The plasmasphere is 
confined to closed magnetic field lines and extends to approximately 4 earth radii at the equator 


2.3 The F2 Peak, Protons, and the Plasmasphere 

Returning to Fig. 3, we see that the majority of the ionospheric plasma peaks at altitudes 
above the FI layer, known as the F2 peak. Although this layer dominates the ionosphere, it 
is not formed primarily by direct photo-ionization occurring at this altitude, but instead is 
the result of other, subsequent processes, including for example, the upwards diffusion of 
0 + ions to higher altitudes where they do not encounter significant numbers of molecules 
with which to charge exchange. In addition, neutral wind drag-driven ion flows are believed 
to be another significant process that contributes to the F2 ionization peak properties at mid- 
and low-latitudes, as discussed by Rishbeth (1967). Finally, upward plasma transport due 
to electric fields, particularly during the daytime, also produces changes in the F2 peak, as 
discussed in Sect. 4. 

Note further in Fig. 2 that the lighter hydrogen ion, H 1 , does not become dominant until 
much higher altitudes. Within the “topside” ionosphere (i.e., above the F2 peak), charge 
exchange between oxygen ions and hydrogen atoms is common, accounting for all of the 
protons of ionospheric origin at mid and low latitudes: 

0+ + H -* H+ + O 

For moderate solar activity, at altitudes above about 1000 km, the H 1 ion dominates, and 
this region is often referred to as the protonosphere. The amount of H 1 present varies con- 
siderably between day and night and is also highly dependent on solar activity. For example, 
Heelis et al. (2009a, 2009b) showed that during extreme solar minimum conditions, the 
crossover altitudes of 0 + and H - ions could be as low as 450 km at night. He + also in- 
creases in importance with altitude, although, except possibly during geomagnetic storms, 
it is typically minor compared to H + . 

The large concentrations of light ions extend to altitudes of about 4 earth radii at the 
equator, as shown in Fig. 4. This high-altitude region is called the plasmasphere and is basi- 
cally a torus that surrounds the mid- and low-latitude portions of the earth and is considered 
an upward extension of the ionosphere at magnetic latitudes below about 60 degrees. The 
plasmaspheric plasma is trapped on closed magnetic field lines that rotate with the Earth 
and thus the plasmasphere co-rotates with them. (The earth’s magnetic field is discussed in 
the next major section.) The outer plasmasphere is composed mostly of hydrogen ions with 
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singly ionized helium accounting for approximately 10-20 % of the plasma (e.g., Craven 
et al. 1997) and has temperatures of roughly 1 eV (or 11,600 K). The plasmasphere ter- 
minates at a sharp decrease in density known as the plasmapause. At lower altitudes, it is 
sometimes associated with a sharp decrease in the topside ionosphere with increasing lati- 
tude that is referred to as the trough. The trough occurs near 60 degrees magnetic latitude, 
and is also used to define the equatorward edge of the high-latitude ionosphere (see review 
by Rodger 2008). 

2.4 Temperature of the Near-Earth Plasma 

In addition to its number density, the near-earth plasma is also characterized by its tempera- 
ture. For the photo-ionization process described above, when an atom absorbs a photon and 
emits an electron, the resulting photoelectron typically has an energy of several eV. This 
energy then dissipates via inelastic collisions with the neutral gas. The electron has a higher 
velocity than that of the ion at the same temperature due to its much smaller mass and the 
conservation of momentum. The resulting temperature distribution of each species depends 
on their collisions with other neutral and ionized particles and thus varies with altitude and 
geographic position. 

Typical average values of electron and ion temperatures vary between 350-3000 K, or 
0.03-0.3 eV, as shown in Fig. 5 for conditions at the magnetic equator at dawn, noon, sun- 
set, and midnight as computed from the IRI model. In the lower ionosphere (<120 km), the 
average values of the electron, ion, and neutral gas temperatures (T e , 7) and T n ) are sim- 
ilar due to the increased number of collisions with the over-abundant neutrals. The values 
of T e , T h and T n then diverge with increasing altitude. The electron temperature is highest 
in the early morning. The neutrals, rotating into sunlight, have yet to warm and the neutral 
atmosphere is still contracted from the colder, nightside temperature. Hence, the electron- 
neutral collisions are less frequent. During the day, the neutral atmosphere expands, and 
the electron-neutral collision frequency thus increases at the higher altitudes, reducing the 


Plasma (T e ,T.) and Neutral (T n ) Temperatures 

IRI-2011 Model (0° lat, 0° Ion) 
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Fig. 5 IRI model calculations showing the plasma and neutral temperatures at the magnetic equator at four 
different local times 
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electron temperature. The electron temperature further decreases during the day as it equili- 
brates with the ion temperature due to electron-ion collisions. The temperatures are largest 
at the higher altitudes where the neutral and plasma densities are low. The enhanced elec- 
tron temperature near 300 km during the day results from photoelectrons produced within 
the upper altitude portion of the layer where the neutral atmosphere is less dense. At night, 
there are no newly-formed photoelectrons and the electron, ion and neutral temperatures 
quickly converge to similar values via collision processes. See Schunk and Nagy (2009) for 
a detailed discussion of the neutral and plasma temperature variations. 

At higher latitudes, precipitating particles, magnetospheric electric fields, and electric 
currents also play significant roles in heating the ionosphere. Frictional heating (Joule heat- 
ing) between the ions and neutrals in particular contributes significantly to ion and neutral 
heating at high latitudes, and plasma instabilities at low altitudes (e.g., two-stream instabili- 
ties, particularly in the auroral zone) may produce elevated electron temperatures locally. 

2.5 Metallic Ions 

Another source of ions to the earth’s near-earth plasma results from the influx of meteoroids. 
The earth is laden by ~10 8 kg/year of meteoritic material (e.g., Ceplecha 1992), depositing 
dust, grains, and atoms released by ablation. As these particles settle due to gravity, metallic 
atoms may undergo charge exchange to form metallic ions such as Fe + , Mg + , and Na + , and 
photo-ionization may also occur (Kopp 1997). 

The main source of metallic ions in the ionosphere, however, is impact ionization by 
meteors striking the atmosphere. This mechanism directly produces ionized gases at the 
meteor head, which occurs when the kinetic energy of the incoming meteoroid is sufficient 
to remove electrons from the neutrals (e.g., Vondrak et al. 2008). Whether a meteoroid will 
undergo impact ionization depends on its mass, velocity, and angle of attack with respect 
to the atmosphere. Such ionization is observed through radar echoes of meteor trails, as 
they interact with the atmospheric gas primarily between 80-110 km, forming cylinders of 
ionizing gas often orders of magnitude denser than that of the surrounding plasma (e.g., 
Dyrud et al. 2001; Close et al. 2002; Janches et al. 2009). An example of such radar “head 
echoes” associated with the ionized gas at the head of an ablating meteor is shown in Fig. 6 
(Chapin and Kudeki 1994). Notice in particular how the head echo extends for over 15 km in 
altitude. The head echoes are followed by a radar trail of plasma irregularities with organized 
scatter that, in this case, lasted for several minutes as they were formed near the magnetic 
equator. Such head echoes are commonly detected by ionospheric radars at all latitudes and 
illustrate the prevalence of meteor-produced ionization at low altitudes (e.g., Mathews et al. 
2008). 

The glowing gas associated with ablation of larger meteoroids produces the familiar 
“shooting stars” that are observed in the night sky at all latitudes and which are most preva- 
lent in the pre-dawn hours when the earth’s proper motion directly sweeps through the in- 
terplanetary medium. For such emissions to be visible to the unaided eye, the mass of the 
meteoroid is generally on the order of milligrams, whereas “radar” meteors have consider- 
ably smaller masses, on the order of 1-100 micrograms, and their ablation optical signature 
is generally not visible to the unaided eye. Ablation of these smaller meteoroids are much 
more frequent and are believed to be a more significant source of metallic ions, at least 
compared to that produced by the larger meteoroids (D. Janches, personal communication, 
2012 ). 

Metallic ions in the earth’s ionosphere have recombination lifetimes for neutralization 
that are very long. Thus, the long lasting metallic ions generally undergo subsequent mo- 
tions after they are created. For example, as the metallic ions settle in the lower portion of 
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Fig. 6 Example of radar backscatter measurements of a meteor ablating in the upper atmosphere. The diag- 
onal line on the left hand side corresponds to the head echo, which creates a plasma density envelope as the 
meteoroid penetrates the atmosphere. The broad scattering region to the right of the diagonal line corresponds 
to plasma irregularities in the meteor trail (Chapin and Kudeki 1994) 


the upper atmosphere, their distributions are transformed by neutral wind shears to form 
“sporadic-E” layers, as discussed in Sect. 4, while others contribute to charged aerosols that 
are responsible for strong radar echoes at lower altitudes in the mesosphere. The long-lived 
metallic ions may also be transported over long distances via upper atmospheric winds and 
electric fields and have been measured by satellites at altitudes as high as 1000 km (Gre- 
bowsky and Aiken 2002), far from the altitudes where ablation is greatest — near and below 
100 km. 

2.6 Daily Variation of the Mid- and Low-Latitude Ionosphere 

Prior to discussing the coupling of the ionosphere to the neutral gas and its dynamics, model 
computations are provided to illustrate how the ionosphere varies, at least in general sense, 
with both local time and altitude. Such a representation is provided in Lig. 7 for both solar 
minimum and solar maximum conditions. These computations are from the 2011 Interna- 
tional Reference Ionosphere (IRI) as discussed in Bilitza et al. (2011), and show conditions 
for all local times along the geomagnetic equator at the equinox. 

The model ionosphere in Lig. 7 shows how the plasma densities are enhanced during 
solar maximum conditions for all local times, with the peak of ionization (the L2 peak) 
moving to significantly higher altitudes compared to solar minimum conditions. Lor both 
solar minimum and solar maximum, the ionosphere is considerably more dense at higher 
altitudes (>250 km) for most of the day and well into the night, exhibiting its lowest plasma 
number densities before dawn. The plasma density at the lower altitudes (<200 km) is very 
low during the night due to recombination losses of the dominant molecular ions. The ap- 
parent decrease in density of the base of the L-region after sunset, which is most pronounced 
during solar maximum, is due to vertical transport and is discussed in Sect. 4. 
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Electron Density Along Magnetic Dip Equator 
IRI-2007 Model 


Solar Minimum (Sept 2008) Solar Maximum (Sept 2000) 

Geo Longitude (deg) Geo Longitude (deg) 
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Fig. 7 IRI model calculations of electron density along the magnetic dip equator for 24 hours of local time 
at solar minimum {left) and solar maximum {right) conditions at the equinox 


Another representation of this same equatorial plasma density distribution is shown in 
Fig. 8, in which the density is shown in a global distribution around the planet along 
the magnetic equator, again computed from the IRI. The panel depicts the ionosphere to 
2000 km and is to scale with the earth’s diameter. Seen as a narrow black circle, a very thin 
ionosphere-earth cavity (i.e., between the lower edge of the ionosphere and the earth’s sur- 
face) is present, whose dimensions increase at night as the lower ionosphere goes away due 
to recombination effects, as discussed above. This cavity is particularly important for the 
excitation of Schumann resonances (low frequency electromagnetic waves that propagate 
around the globe) and other phenomena, as discussed by Simoes et al. (this issue). 

To illustrate the diurnal variations of the ionosphere at a single location, the measured 
plasma density from 50-500 km is plotted for four consecutive days in Fig. 9 as observed 
by the Arecibo observatory in Puerto Rico (Figure courtesy N. Aponte). The day numbers at 
the bottom of the plot correspond to local noon. Clearly, the data show that the lower altitude 
(<250 km) ionosphere turns on very quickly after sunrise and then it recombines abruptly at 
sunset. During the night, the ionospheric density is hence substantially reduced at altitudes 
up to about 250 km, where the lower edge of the F-region can be very “ragged”, influenced 
perhaps by transport effects. During the day, the F-region plasma density is fairly uniform, 
extending into nighttime where it begins to decrease at the higher altitudes after midnight. 
The very faint, thin lines visible at night near 110 km are due to sporadic -E layers and will 
be discussed further on below. 

Finally, we show IRI model computations illustrating how the plasma temperature varies 
with local time and altitude in Fig. 10. Note how the electron temperature increases dramat- 
ically in the early morning, as discussed above. 


3 Coupling of the Near-Earth Plasma and the Neutral Atmosphere Gases 

In this section, we first describe how the ionized gases are organized by the ambient mag- 
netic field and then discuss how the ionospheric charged constituents are coupled to the neu- 
tral gas via collisions, which vary markedly with altitude. We describe how this coupling is 
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Fig. 8 IRI model calculations of the plasma density at the magnetic equator, presented in the altitude region 
scaled to the earth’s radius 
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Fig. 9 Incoherent scatter measurements of the plasma density for four consecutive days at the Arecibo 
Observatory in Puerto Rico (Figure courtesy of N. Aponte) 
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Fig. 10 IRI model calculations 
of the electron and ion 
temperatures along the magnetic 
equator for a 24 hour period 


Plasma Temperatures Along Magnetic Dip Equator 

IRI-2007 Model, 22 Sep 2004 
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Fig. 11 The Earth’s dipole 
magnetic field. The field lines are 
nearly vertical at very high 
latitudes yet nearly horizontal at 
low latitudes. The magnetic field 
line geometry has profound 
effects on the near-earth plasma. 
The magnetic poles are offset 
from the geographic poles 



represented by the mobility and conductivity tensors, which form the basis of a large number 
of physical processes that govern the behavior of the partially ionized gas environment. 

3.1 The Earth’s Magnetic Field 

Unlike the neutral atmospheric gases, the ions and electrons that comprise the near-earth 
plasma are compelled by Lorentz forces to gyrate about the earth’s magnetic field lines. 
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Earth's Magnetic Equator 



Longitude 


Fig. 12 The lower panel shows the location of the earth’s magnetic equator. The upper panel shows the 
model magnetic held strength along the magnetic equator at 100 km and 500 km altitude 


Thus, an appreciation of the earth’s magnetic field is essential to understand the behavior of 
the near-earth plasma. 

The Earth’s magnetic field is essentially that of a dipole whose axis is tilted with respect 
to the earth’s axis of rotation by approximately 10 degrees, as shown in Fig. 11. Conse- 
quently, the magnetic equator is also offset from the geographic equator. This figure shows 
how the main field lines are nearly vertical at high latitudes, stretching into the deep magne- 
tosphere, yet remain closer to the earth, subtending much smaller angles with respect to the 
horizontal below about 45 degrees latitude. Indeed, at latitudes of roughly 60 degrees and be- 
low, the field lines “close”, meaning that they extend to the conjugate hemisphere at roughly 
the same magnetic latitude. This condition strongly influences the near-earth plasma, as we 
shall show. 

The magnetic field strength at 100 km altitude near the poles is ~0.65 G, and decreases 
to about half that value at the magnetic equator, which is defined by the location on the earth 
where the magnetic field vector is strictly horizontal. In addition to being tilted, the dipole 
is slightly distorted and hence the magnetic field strength near 100 km altitude may be as 
low as 0.25 G in the Peruvian sector but near 0.4 G in the Indian sector at the magnetic 
equator. These differences strongly influence currents and longitudinal effects of plasma 
flow. Figure 12 shows the geographic location of the model magnetic equator, as well as 
the strength of the magnetic field at the equator. The declination of the magnetic field at 
the equator is generally perpendicular to the magnetic equatorial plane. This declination 
introduces longitudinal influences on plasma drifts and can also be important for magnetic 
flux tubes whose “footpoints” in the northern and southern hemispheres may be significantly 
displaced in longitude. 

As we shall show, the earth’s intrinsic magnetic field has a profound effect on the plasma 
with respect to directing E x B drifts and currents, as well as connecting the near-earth 
plasma with energy and momentum sources in the distant magnetosphere. 

3.2 Coupling of the Ionized and Atmospheric Gases: Collisions and the Mobility and 
Conductivity Tensors 

Overview. We demonstrated in the previous section how the ionized and atmospheric gases 
co-exist in the same volume particularly in the altitude range of 90-1000 km. We now dis- 
cuss how the ionized and neutral gases are coupled, which sets the stage for later sections 
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Fig. 13 Altitude profiles of the 
electron-neutral collision 
frequency and the ion-neutral 
collision frequency. Profiles of 
the mean ion gyrofrequency and 
electron gyrofrequency are also 
shown 


Wallops Island, Virginia, USA 
Midday 



where we show how the partially ionized fluid responds to different forcing functions, pri- 
marily those due to electric fields and winds. In particular, we demonstrate that, although the 
ionized component is controlled by the earth’s magnetic field, both the ionized and neutral 
gas motions are subject to collisional interactions between the ions and neutrals that vary 
with altitude. These interactions are specified by the mobility and conductivity tensors, as 
shown below. 

Gyro and Collision Frequencies. Charged particles in the ionosphere gyrate about the 
ambient magnetic field lines at their cyclotron frequencies, Qj, given by 

£lj=e\B\/mj (1) 

where e is the electric charge, |B| is the magnetic field strength, and m ; is the mass of 
the species which is designated by j . For reference, the gyro frequency for 0 + ions in a 
magnetic field of 0.5 Gauss is 48 Hz, whereas the gyro frequency for electrons is 1.4 MHz 
for the same magnetic field strength. 

Since the upper atmospheric gases are dense (compared to the plasma) and only weakly 
ionized, collisions between ion and neutrals significantly influence the individual and bulk 
motions of the plasma. Because of these collisional effects, the statistical collective behavior 
of the particles may be defined using the equations of hydrodynamics. Various expressions 
have been put forward for the ion-neutral collision frequency, v in , and the electron-neutral 
collision frequency, v en , that depend on the neutral density, neutral and ion masses, and the 
temperatures (e.g., Banks and Kockarts 1973). 

For the lower altitude regions of the ionosphere, model ion-neutral and electron-neutral 
collision frequency profiles are shown in Fig. 13 for a mid-latitude location at mid-day. It 
is clear that below 130 km, the ions are collision dominated by the atmosphere whereas the 
electrons are able to execute their gyro motions about the magnetic fields and consequently 
undergo E x B drifts across magnetic field lines. The lower ionosphere, where v in / Q, > 1 
and v en / Ff, <$; 1 , constitutes the dynamo region where strong currents flow, as discussed in 
the next section. 
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Mobility and Conductivity Tensors. Because the plasma distribution function in the pres- 
ence of a magnetic field is highly anisotropic, we define a mobility tensor, JZ, which relates 
the velocity (V) and force (F) vectors: 


where 


V = M-F 


M = 


Mi 

-M2 

0 


M 2 

Mi 

0 


0 

0 

Mo 


( 2 ) 


(3) 


This tensor has been defined in a geometry such that the magnetic field, B, is parallel to 
the z-axis. For single particle motion subject to an applied force, we have the following 
mobilities: 


fjL 0 j = l/mvj n (parallel mobility) 

Mi/ = (1/ mvj n )(vj n ) 2 / (Q, 2 + v 2 n ) (Pedersen mobility) 

M2 j = (l/mv jn )(v j „Qi)/(Q, 2 j + vj fl ) (Hall mobility) 

The profiles of these three mobility components are plotted in Fig. 14 for three different 
locations, one at the magnetic equator (Jicamarca, Peru), one at high latitudes (Poker Flat, 
USA), and one at mid-latitudes (Wallops Island, USA). The mobilities shown here are com- 
bined for ions and electrons, and use the mean molecular weight for the ions. Note that the 
Pedersen mobility peaks at around 125 km and the Hall mobility peaks near 105 km. The 
mobilities are higher at the lower latitudes due to the weaker magnetic field strength, since 
they are roughly proportional to ~1/ Qf. 

We now define a conductivity tensor, a, that relates the electric field and current density: 

J = W E (4) 

which can also be expressed in terms of the difference of the ion and electron motions: 

J = Ne(\i — \ e ), (5) 

where N is the plasma density and e is the electric charge. 

The conductivity may be expressed in terms of the mobilities such that a = 
Ne 2 (ii i + Me)- The components of the conductivity are thus: 

a Q = Ne 2 (/i oe + Mo;) (parallel conductivity) 
o\ = Ne 2 (/i ie + Mu) (Pedersen conductivity) 
er 2 = /Vc 2 (m 2 <- — M2i) (Hall conductivity) 

The ionospheric conductivity depends linearly on the ambient plasma density and the cou- 
pled neutral and ionized gases via the anisotropic mobilities as organized by the magnetic 
field. Equation (4) may now be expressed more fully: 

J = cr | E pcrp + (To b' x E perp + er 0 E par (6) 

where b' is the unit vector along the magnetic field direction and E perp is the perpendicular 
electric field in the earth-fixed frame, such that E perp = E perp + U x B, where U is the neutral 
wind. This expression shows that if there is an electric field perpendicular to the magnetic 
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Fig. 14 Altitude profiles of the Pedersen and Hall mobilities at Wallops Island, Virginia; Poker Flat, Alaska; 
and Jicamarca, Peru 


field, B, the current will have two components — one along the electric field but perpendic- 
ular to B (Pedersen) and one perpendicular to both the electric field and to B (Hall). Their 
magnitudes differ by v, / £2, and the angle, a, between them is given by a = tan -1 (£2,/u,). 

Model calculations of the Hall and Pedersen conductivity in the E-region at different 
latitudes are presented in Fig. 15. (The parallel conductivity for all locations is orders of 
magnitude higher than that of the corresponding Hall and Pedersen conductivity and is not 
shown in the figure.) The top two panels show the plasma density at noon and the Hall and 
Pedersen conductivities calculated using the mobilities as shown in Fig. 14. The lower panels 
show the plasma density and conductivity at 2100 FT for Poker Flat, Alaska for both quiet 
(left) and active (right) conditions. The active auroral conditions were selected to provide 
plasma densities of 10 5 /cm 3 at 1 10 km, corresponding to a moderate auroral arc. Clearly, for 
conditions where there is no aurora, the high latitude conductivity at night is very low. The 
thermal plasma created by the auroral particles provides a very strong conductivity which is 
important for current closure, Joule heating, and other important high latitude effects, as we 
will discuss later on. 


4 The Dynamics and Morphology of the Ionosphere at Low and Middle Latitudes 

In Sect. 2, we showed how the earth’s ionosphere at low and middle latitudes is formed 
via solar photo-ionization and subsequently varies with altitude and local time, revealing 
important differences between the daytime and nighttime plasma number density and tem- 
perature. In this section, we show how this ionization is significantly altered by dynamics. 
The dynamics are driven primarily by neutral forcing (winds) of the ions that subsequently 
set up internal polarization electric fields. In the next section, we will show how the opposite 
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Fig. 15 Altitude profiles of the Hall and Pedersen conductivities for representative plasma density profiles 
at Wallops Island, Virginia, USA ( upper left), and Jicamarca, Peru ( upper right). The lower panels show the 
Hall and Pedersen conductivities at Poker Flat, Alaska, USA, for quiet conditions {left) and moderate auroral 
conditions (right) 


occurs at high latitudes, i.e., electric fields drive E x B plasma motion via large amplitude, 
externally imposed, magnetospheric fields that, through ion-neutral collisions, subsequently 
drive neutral winds. 

4.1 Dynamo Currents and the Equatorial Electrojet 

As the solar EUV heats the neutral upper atmosphere, the neutral gases generally expand 
away from the sub-solar point, creating a system of upper atmosphere tidal motions or winds. 
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Fig. 16 Current patterns set up by the atmospheric dynamo for daytime {left) and nighttime {right) hemi- 
spheres (after Matsushita 1965) 


The neutral motions drag ions across magnetic field lines via ion-neutral collisions, creating 
a very small charge displacement with the electrons that subsequently sets up relatively 
weak, ionospheric electric fields. At the lower altitudes (<130 km), the ion-neutral collision 
frequencies become so large that the ions are unable to complete their gyro orbits whereas 
the electron cyclotron orbits, due to their much smaller mass, are not restricted and hence, 
remain “magnetized”. Due to the differential drifts of the ions and electrons at these low 
altitudes, a system of global scale, horizontal current patterns develops, in which internal, 
polarization electric fields are generated to restrict the current flow to be non-divergent. This 
worldwide system of currents is called the “solar quiet” (or Sq) dynamo which form current 
loops at mid-latitudes which close at the magnetic equator, as shown in Fig. 16 (Matsushita 
1965). 

The strongest current driven by the atmospheric dynamo is the equatorial electrojet. This 
current system exists as a result of both the dynamo action at all latitudes and the enhanced 
conductivity set up by the horizontal magnetic field geometry at the magnetic dip equator. At 
this location, in the altitude range of roughly 95-1 15 km, significant zonal (east-west) elec- 
tron drifts are set up in the following way. Vertical E x B drifts due to the ambient eastward 
electric fields (driven by the solar forcing discussed above) combined with the north-south 
magnetic field are inhibited for ions because the ion-neutral collisions are sufficiently large 
such that Vi = 0. Hence, only the electrons are magnetized and the resulting vertical charge 
separation between the electrons and ions establishes a vertical polarization field, E p . This 
polarization electric field, in turn, creates a horizontal E p x B drift which reinforces the 
small, existing dynamo current and thus creates the strong electrojet current. This process 
is depicted in Fig. 17. The equatorial electrojet is strongest within a degree of the magnetic 
equator, where the magnetic field vector is precisely horizontal, thus enabling the vertical, 
polarization electric field described above to be set up. This current can also be explained 
in terms of the Cowling conductivity, which is an enhancement, by a factor of er| /cr x , of the 
Pedersen conductivity at the magnetic equator (e.g., Kelley 2009). 

At night, the zonal electric field reverses direction and thus the polarization electric field 
is downward. The resulting E x B drift again enhances the electrojet current since the current 
is also of opposite direction at night. The magnitude of the nighttime current densities is 
much lower than during the daytime due to the greatly reduced ambient electron density. 

An example of measurements of the polarization electric field, current, and plasma den- 
sity measured on a sounding rocket launched in the daytime electrojet from Brazil is shown 
in Fig. 18a. In this figure, the resulting electrojet plasma drifts are calculated independently 
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Fig. 17 The zonal electric fields create the vertical polarization fields (E p ) at the magnetic equator in day- 
time {left) and nighttime (right). The E p x B drifts drive the daytime and nighttime electrojets in opposite 
directions in accordance with the sign of E p 


using both the measured electric fields (i.e., the E x B drift) and the measured current and 
plasma density (i.e., \ e = J/qN e , assuming V, = 0), which show comparable results (Pfaff 
et al. 1997a). The discrepancy between these two drift profiles may be due to the neutral 
wind contribution, which was not measured in this experiment and hence not taken into 
account. 

When the relative electron-ion drift speed exceeds approximately the acoustic veloc- 
ity, collisional two-stream waves are excited (Farley 1963; Buneman 1963). These intense, 
meter-scale waves are observed in electric field wave measurements recorded by sounding 
rockets, for example in the equatorial electrojet off the coasts of Peru and Brazil (Pfaff et al. 
1987, 1997b). The altitude- spectrogram in Fig. 18b shows an example of such two-stream 
waves, represented by the broadband emissions extending to higher frequencies between 
102 and 110 km, where the instability threshold conditions are met. The intense, lower fre- 
quency waves observed in particular at the lower altitudes are driven by the gradient-drift 
instability and are set up by the vertical (positive upwards) DC electric field in concert with 
the upward ambient plasma gradient. The lower edge of the daytime ionosphere is unstable 
at long wavelengths in the electrojet region, as evidenced by structure in the plasma den- 
sity profile at 105 km altitude and below, shown in the figure and discussed in Pfaff et al. 
(1987). Indeed, the lower ionosphere in both the daytime and nighttime electrojet regions is 
typically highly unstable to these plasma instabilities (e.g., Fejer et al. 1975). Distinct, km- 
scale waves dominate the spectrum when the gradient-drift instability is active, propagating 
zonally along the electrojet direction. These large scale waves have been shown to change 
directions between daytime and nighttime conditions as the ambient polarization electric 
field changes sign (Kudeki et al. 1982, 1987). 

4.2 The Fountain Effect: The F-region Dynamo, and Plasma Drift Patterns at Low 
Latitudes 

Given the formation of the ionosphere by photo-ionization from solar photons and the con- 
trolling role of the solar zenith angle in the resulting plasma production, it would be logical 
to expect the ionosphere to have maximum number density at the sub- solar point and that 
the variations in this density would be fairly predictable with local time, latitude, and season. 
In fact, this is not the case. The daytime ionospheric density does not peak at the sub-solar 
point, but actually has two distinct maxima about ±15 degrees from the magnetic equator 
which can be quite pronounced. 


Springer 


The Near-Earth Plasma Environment 


43 


NASA 21.111 (Downleg) -- Alcantara, Brazil 

21 September 1994 10:56:00 L.T. 

Zonal DC E Field Vertical DC E Field Current Density Plasma Density Electron Drift Velocity 




(a) 

Flight 33.027 -- Punta Lobos, Peru 
March 12, 1983 -10:35 L.T. 
Upleg 



Fig. 18 (a) Sounding rocket observations of the altitude profile of the DC electric field, current density 
and plasma density in the daytime equatorial electrojet at Alcantara, Brazil. The panel on the right shows 
independent calculations of the electron drift velocity using (solid) the current density and plasma density 
and ( dotted ) the polarization electric field (Pfaff et al. 1997b). (b) Sounding rocket observations of electric 
field waves displayed in an altitude sonogram ( right panel) in the daytime equatorial electrojet at Punta Lobos, 
Peru. The color levels correspond to over 40 dB in wave power, with red denoting stronger amplitudes than 
blue. Simultaneous measurements of the unstable plasma density are shown in the left hand panel (Pfaff et al. 
1987) 

These ionization peaks and the equatorial valley in between them were initially observed 
in ionosonde data and eventually explained as a consequence of the eastward electric field 
in the daytime, low latitude ionosphere (Appleton 1946; Martyn 1947). As a result of this 
electric field, the E x B drift forces plasma upwards at the magnetic equator which then, 
via diffusion and gravity, “slides” down the magnetic field lines to form the two maxima 
away from the magnetic equator, peaking at ~300-400 km altitude, as illustrated in Fig. 19. 
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Electron Density, SAMI2 Model 
22 Sep 2009, SLT = 15 hr, Long. = 205° E 
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Fig. 19 SAMI2 model calculations versus latitude and altitude of the plasma density for 1500 SLT at 205° 
East. The upward E x B drift at the magnetic equator is driven by the eastward electric field, and there is 
subsequent flow downward along the magnetic field lines 


Fig. 20 SAMI2 model 
calculations of the global total 
electron content (TEC) at 
22:32 UT, showing the equatorial 
anomaly enhancement in the late 
afternoon. The center of the 
globe corresponds to 17:44 LT 
(Figure courtesy J. Fluba) 



Referred to as the so-called “fountain effect” or the equatorial ionization anomaly (EIA), 
this phenomenon represents a significant, fundamental mechanism that characterizes the 
near-earth plasma environment. Model total electron content calculations for moderate solar 
conditions are shown in Fig. 20 that were calculated with the SAMI2 model (Huba et al. 
2000) and show how the anomaly “wings” dominate the late afternoon ionosphere at low 
latitudes. 

Another process at work at low latitudes that drives the plasma is the so-called “F-region 
dynamo”. At high altitudes in the low latitude F-region during the daytime, wind-driven 
field aligned currents flow along the magnetic field lines between the F and E-regions. At 
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night, however, the E-region ionosphere is devoid of plasma due to recombination effects 
discussed previously, and hence, the “feet” of the magnetic field lines, whose apex altitudes 
are in the F-region, are no longer in high conductivity E-region. As a result, large-scale 
electric fields are set up in the F-region which drive a low-latitude dynamo in the F-region 
at night in order to maintain current continuity (e.g., Rishbeth 1971, 1997; Heelis 2004). 

An enhanced vertical drift in the F-region after sunset is a characteristic of the low lat- 
itude ionosphere, whose intensity exhibits seasonal, longitudinal, and solar cycle depen- 
dencies, tending to be largest near the equinoxes at solar maximum (Scherliess and Fejer 
1999). This “pre-reversal enhancement” has been attributed to processes also related to dif- 
ferences in the conductivity between the E- and F-regions near the terminator and is gen- 
erally attributed to a polarization electric field set up to maintain current continuity near 
the dusk meridian (e.g., Heelis et al. 1974; Farley et al. 1986; Haerendel and Eccles 1992; 
Eccles 1998). 

Although the winds and polarization electric fields can be quite variable day-to-day, an 
average statistical pattern emerges from in situ and radar measurements of meridional (ra- 
dially outward) drifts at low latitudes which is used in a number of models to represent 
the low latitude electrodynamics. Measurements of the meridional drifts near the magnetic 
equator from the ROCSAT satellite (Fejer et al. 2008) averaged for moderate solar activity 
conditions are shown for all longitudes and local times in Fig. 21. Note how the diurnal 
variations have a marked dependence on longitude, with the strongest enhancements after 
sunset occurring in the longitude sector of —120 to +60 degrees East at the equinox. 

Because the mid- and low-latitude ionosphere is a “closed” system (i.e., the magnetic 
field lines that encompass this region are all closed), it may appear straightforward to de- 
velop a self-consistent electrodynamic model based on the wind driven E- and F-region 
dynamos, even when electric field mapping and variations of magnetic field declination 
with longitude are taken into account. However, the dynamo-driven drifts are not the whole 
story. Tidal effects (discussed next) propagating from below also have profound effects in 
the ionosphere, and account for the 4-cell pattern in longitude in the morning vertical drifts 
shown in Fig. 21. In addition, neutral forcing from high latitudes as well as penetration elec- 
tric fields of magnetospheric origin at all latitudes due to strong magnetic activity will also 
affect the global drift patterns. Finally, upwards propagating electric field contributions from 
tropospheric storms are also believed to contribute to the drifts at any given time, particu- 
larly at night when the plasma density is greatly reduced and the thunderstorm-generated 
electric fields couple more efficiently to the ionosphere. 

4.3 Effects of Tides on the Plasma Density and Electric Fields 

A very important contribution to the low and mid-latitude upper atmosphere is neutral forc- 
ing from below. The strongest and most regular forcing is due to solar thermal tides which 
exist over a wide range of harmonics related to the earth’s 24 hour rotation period. Given 
the strong effect of neutral winds on the low and mid-latitude electrodynamics, it is not 
surprising that tidal forcing might have profound effects on the near-earth plasma. The sig- 
nificance of these tidal effects and the fact that they are well-organized, for example with 
longitude, have received renewed interest in the space physics community as new data sets 
have recently become available. 

Solar tidal effects have been observed in the ionosphere in both the electrodynamics and 
the plasma density. The actual mechanisms by which the tidal signatures are manifest are 
still being debated. Apart from migrating tides which are responsible for creating day-night 
variations in the ionosphere, the most distinct tidal signature is attributed to the diurnal east- 
ward wavenumber-3 (DE3) non-migrating tide which produces very clear longitude effects 


Springer 


46 


R.F. Pfaff 


Fig. 21 Average vertical drift 
velocities recorded by the ion 
drift meter on the ROCSAT 
satellite, as a function of local 
time, longitude, and season, for 
quiet times during moderate solar 
flux conditions (Fejer et al. 2008) 
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in the upper atmosphere. An example of how this tidal structure effects the low latitude, 
E-region electric fields is shown by England et al. (2006), who compared the equatorial 
electrojet strength as measured by the CHAMP, SAC-C, and Oersted satellite magnetome- 
ters and demonstrated that they each display the same 4-peak longitudinal structures. An 
example of enhancements in the total electron content at evenly-spaced longitudes at low 
latitudes for the 12-14 LT sector is shown in Fig. 22 using observations from the COSMIC 
satellites. (Figure provided courtesy of C. Lin. See Lin et al. 2007, for details.) These au- 
thors argued that the tides influence the F-region plasma by affecting the dynamo electric 
fields which subsequently modulate the vertical plasma transport. To demonstrate this, these 
authors used the same COSMIC data to create cross section plots of the equatorial ionization 
anomaly for a peak and a valley period to show that the EIA extended to higher altitudes and 
higher latitudes in periods of peak tidal forcing, a behavior which can be clearly linked to 
larger driving electric fields. That the equatorial electric fields show this same wave-4 tidal 
signature is also evident in the 4-peak longitude variations of the vertical plasma drifts evi- 
dent in the morning (9-1 1 LT) period in Fig. 21. A review of the effects of the non-migrating 
tides on the near-earth plasma is provided by England (2011, this issue). 

Airglow enhancements in the crests of the equatorial anomaly peaks at 20 LT provides 
another representation of this non-migrating D3 tide phenomenon, as shown in Fig. 23 (Im- 
mel et al. 2006). This figure shows a reconstruction of 30 days of observations gathered 
with the FUV imager on the IMAGE satellite. The southern “band” is introduced as a mi rror 
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Fig. 22 Average TEC observations recorded by sensors on the COSMIC-2 satellites for the 1200-1400 LT 
sector as a function of latitude and longitude ( upper panel). The lower panels show cross sections versus 
altitude and latitude corresponding to a minimum {left) and peak {right) TEC region, illustrating how the 
equatorial anomaly responded to different electric fields in each case (figure courtesy C. Lin. See also Lin 
et al. 2007) 



Fig. 23 Reconstruction of nighttime ionospheric emissions for March 20-April 20, 2002, as measured with 
the IMAGE satellite (Immel et al. 2006), which shows the average location and brightness of the equatorial 
anomaly. (Because of poor sampling of the southern hemisphere, its emissions are represented by a mirror 
image of the northern anomaly.) The image represents data at 2000 LT and includes an overlay of white 
dashed contours which represent model temperatures at 115 km 


image of the upper, northern band. The white dashed contours shown in this presentation 
represent model temperature variations due to the upward propagating tidal oscillations, 
which correlate nicely with the observed peaks in the equatorial anomaly. 
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Fig. 24 Backscatter radar power 
at 50 MHz as recorded on the 
night of March 6, 1983. Darker 
shades correspond to more 
power, which is depicted on a 
logarithmic scale (LaBelle 1985; 
Kelley et al. 1986) 
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4.4 Large Scale Plasma Depletions at Low Latitudes (“Spread-F”) 

At low latitudes at night, major distortions of the near-earth, ionospheric plasma may occur, 
often of planetary-scale proportions. This phenomenon may be characterized by a series of 
depleted magnetic flux tubes or wedges that stretch between magnetic conjugate locations in 
the northern and southern hemispheres. Sometime referred to as “bubbles”, these elongated 
plasma depletions typically have internal, eastward electric fields of 5-20 mV/m (although 
sometimes much greater) and hence undergo significant upward E x B motions at speeds 
much greater than those of the background, ambient plasma (e.g., Aggson et al. 1992). 
The seat of the depleted flux tubes is well established to be the base or lower ledge of 
the nighttime F-region plasma (e.g., near 250 or 300 km). When fully developed, however, 
these depletions can rise to altitudes of over 1000 km and extend significantly in latitude. 
In some cases, primarily off the magnetic equator, narrow plasma enhancements, instead of 
depletions, are observed (e.g., Le et al. 2003; Klenzing et al. 201 1). 

A key feature of the equatorial plasma depletions is that they are themselves unstable, and 
readily display strong electrostatic plasma irregularities with scale sizes that range from over 
100 km to less than one meter. Certain scales contribute to Fresnel and Bragg scatter at high 
frequencies and cause radio scintillations as well as radar and communication outages (e.g., 
see de La Beaujardiere 2004, and references therein). In fact, this general phenomenon was 
first identified by its associated irregularities, which create a spread in the F-region scatter on 
ionosonde traces and is generally referred to as “spread-F” in the literature. Recent reviews 
of large-scale plasma depletions/spread-F may be found in Woodman (2009), Kelley et al. 
(2011), and Makela and Otsuka (2011, this issue). 

An example of equatorial spread-F is shown in the Jicamarca backscatter radar map for 
3-m irregularities in Fig. 24 (Kelley et al. 1986; LaBelle 1985). Such radar “range-time- 
intensity” plots are similar to observations gathered with a slit camera: the radar beam has 
a fixed, one-dimensional (vertical) orientation through which the two-dimensional, dynam- 
ical ionospheric structures are observed. Thus, the reader must use caution in interpreting 
the image as a two-dimensional map, particularly because the features are not stationary, but 
evolve in space and time. With that precaution, we can nevertheless draw some important 
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Fig. 25 Spatial relationship of 
the backscatter power measured 
during equatorial spread-F by the 
Altair radar at Kwajalein Atoll 
(lower panel) with the density 
measurement by the AE-E 
satellite on August 13, 1978 
(after Tsunoda et al. 1982). The 
satellite trajectory is shown in red 


AE Satellite - Plasma Density Measurements 
Altair Radar (Kwajalein) - Backscatter Echoes 
13 August 1978 


10:41 10:42 10:43 UT 



MAGNETIC EAST DISTANCE FROM ALTAIR - km 


conclusions regarding large-scale ionospheric features and dynamics from this presenta- 
tion. For example, the strong radar backscatter near 400 km at 20:25 UT almost certainly 
corresponds to the F-region ledge which then rises near 20:45 UT, presumably due to the 
pre-reversal enhancement of the eastward electric field, which causes the F-region to rise 
after sunset at low latitudes. The altitude of the ledge then decreases and settles at a range of 
about 300 km during the latter portion of the event shown (i.e., near 21:35 UT). This strong 
scatter that traces out the F-region ledge is known as “bottomside” spread-F. 

The salient feature of Fig. 24 is the spread-F “eruption” shown in the center of the panel 
near 20:45 UT and extends to much higher altitudes, forming so-called “plumes”. In fact, 
in this example, there appear to be 3 or perhaps 4 evenly spaced plumes that emerge from 
the center and backside of the F-region ledge (between 20:45 UT and 21:15 UT) whose 
irregularities constitute “topside spread-F”. Furthermore, note that above about 700 km, 
there are very high altitude structures that display a more diffuse character. Such scatter 
is believed related to smaller scale instabilities that are driven at high altitudes where the 
collision frequencies are negligible. Finally, we note that the narrow, darkened traces near 
100 km correspond to nighttime electrojet irregularities, for which the decrease in intensity 
near 21:00 likely corresponds to the change in the direction of the ambient zonal DC electric 
field, as shown for other examples in Fejer et al. (1975). 

We now consider another radar backscatter example, in this case from the steerable Altair 
radar at Kwajalein, which the AE-E satellite traversed with a plasma density probe, as shown 
in Fig. 25 (Tsunoda et al. 1982). In this case, as shown in the upper panel, the probe clearly 


<£) Springer 



50 


R.F. Pfaff 


100 

10 


1 

400 


S 200- 


-200 

8 

E 4 
> 

E 0 

-4 

9 


> 

E 


0 
-3 
10 5 
10 4 i 
10 
10*-, 
10 1 n 


UTthhimm) 01:38 

ALtfkm) 

ILATCdeq] 


ALtjkm) 452.8 
Tfhrj 9 22 ' 4 


C/NOFS Orbit 688 - June 02, 2008 
VEFI Observations 







_jsa 





MLTli 


1.6 


01:40 

472.7 

21.0 

2.1 


01:42 

495.3 

19.8 

2.6 


Efx § Drift 




Emeridional 




Ion Density 


■ -30 
-40 T 
-50 

O 

-60 o 
Merid (Out) 


Zonal (East) 
Out 


In 

East 


West 


01:44 

520.1 

19.0 

3.1 


Fig. 26 Electric field and plasma density measurements recorded by the C/NOFS satellite during strong 
spread-F density depletions, as shown in the bottom panel. The two panels above the density show the DC 
electric field components and the panel above that shows the corresponding E x B velocity. Notice how the 
zonal component of the electric field, which corresponds to the outward plasma drift component, is anti-cor- 
related with the log of the plasma density. The uppermost panel shows ELF electric field spectra that reveal 
strong irregularities coincident with the density depletions 


observes profound plasma density depletions, decreasing by over two orders of magnitude, 
which correspond precisely to the regions of strong radar backscatter, as shown in the lower 
panel. This example and many others have demonstrated conclusively that strong backscatter 
echoes originate in regions of well-defined, abrupt plasma depletions. 

An example of C/NOFS satellite probe observations of plasma density depletions and 
associated electric fields within the depletions is shown Fig. 26. The data show a series 
of depletions that extend to very low plasma densities and which have associated upwards 
and westward E x B drifts, as measured by the vector electric field detector on that satellite. 
Note how the zonal electric field and the log plasma density are anti-correlated, even at small 
scales, indicative of a quasi-linear relation between these quantities, as shown by Aggson 
et al. (1992) for similar observations using data from the San Marco satellite. The C/NOFS 
data in the uppermost panel of Fig. 26 show electric field irregularities with the associated 
depletions. These are the irregularities that disrupt radiowaves, as discussed above. 

Although the definitive source of the onset and triggering of such low latitude plasma 
depletions is still debated, the rising, depleted flux tubes are believed a consequence of an 
interchange or Rayleigh-Taylor instability (e.g., Woodman and La Hoz 1976). In this case, 
the plasma density at low altitudes (<250 km) is greatly reduced at night so that there exists a 
“heavy” fluid (i.e., the F-region ionosphere) suspended over a “lighter” fluid. The instability 
is particularly susceptible near the magnetic equator because the magnetic field lines in the 
F-region are nearly horizontal and essentially sustain the higher density plasma at higher 
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Jicamarca Combined Incoherent and Coherent Scatter 
24 February 1996 
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Fig. 27 E x B drifts ( white arrows) and Jicamarca backscattered power map for the early evening hours of 
February 24, 1996 (Kudeki and Bhattacharyya 1999). The scattered power corresponds to incoherent scatter 
(brown and green shades), and shows how the lower edge of the ionosphere rises after sunset, which occurred 
at 18:32 LT. The stronger scatter ( red and purple shades) correspond to coherent scatter of the spread-F and 
E-region irregularities. Notice the vortex pattern in the E x B arrows prior to the onset of the high altitude 
plume that begins near 20:30 LT 


altitudes. Small perturbations on the F-region ledge would produce an eastward electric 
field which grows in concert with an associated density depletion, starting the instability 
process. Another factor believed important for irregularity development is that at sunset, 
large pre-reversal, zonal electric fields cause the F-region to rise suddenly to altitudes more 
susceptible to shorter scale instability growth (e.g., Kelley 2009). 

Despite the apparent success of the Rayleigh- Taylor instability explanation for the evo- 
lution of low latitude plasma depletions, it is widely acknowledged that this instability is 
too slow to produce the observed prompt depletions after sunset and an additional trigger- 
ing mechanism is needed to explain their rapid appearance after sunset. Furthermore, the 
presence of a heavy fluid over a light fluid at the magnetic equator occurs every night, so 
some additional instigator is required beyond statistical fluctuations. Recent theoretical anal- 
ysis has shown that spread-F depletions may be triggered whenever eastward neutral winds 
propagate at different speeds than that of the ambient plasma, as discussed by Kudeki et al. 
(2007). These authors show how a vortex in the ambient plasma at sunset (in the plane per- 
pendicular to the magnetic field) contributes to the fact that the plasma and neutral gases 
subsequently move at different velocities. An example of the vortex motion in the plasma 
drift detected by the Jicamarca radar is shown in Fig. 27 (Kudeki and Bhattacharyya 1999). 
Shear instabilities also have been suggested to create the large-scale undulations that ulti- 
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Fig. 28 Simultaneous distributions of total electron content showing gravity waves/TIDS at 23:20 LT on 
May 22, 1998 using GPS data (a) and 630 nm airglow data (b) (Saito et al. 2001) 


mately seed the Rayleigh-Taylor growth of spread-F (e.g., Hysell and Kudeki 2004). This 
shear in the zonal plasma drift at the base of the F-region was observed in altitude profiles 
of the vertical electric field measured on sounding rockets flown from Kwajalein Atoll, as 
reported by Hysell et al. (2006). 

4.5 Acoustic Gravity Waves and Traveling Ionospheric Disturbances 

Acoustic gravity waves have been known for some time to affect the ionosphere and, at some 
level, such oscillations are expected to be ubiquitous in the ionosphere. In general, apart 
from singular events driven by localized phenomena such as earthquakes, two sources of 
gravity waves in the mid and low-latitude ionosphere are customarily studied: (1) upwards- 
propagating gravity waves that originate in the lower atmosphere, such as those gener- 
ated by thunderstorms and other tropospheric energy sources (e.g., Walterscheid 2001), and 
(2) horizontally-propagating gravity waves generated at high latitudes, for example by au- 
roral processes, that subsequently propagate equatorward (e.g., Crowley 1991). The neutral 
motions associated with the gravity waves set the ionosphere in motion where they are ob- 
served as traveling ionospheric disturbances, or TIDs. In a situation similar to that of tidal 
forcing, recent advances in instrumentation and observations have enabled a more in-depth 
analysis of the role on gravity waves/TIDs in near-earth region of geospace. 

Traveling ionospheric disturbances are frequently observed in airglow, radar, and GPS 
measurements. Recently, arrays of GPS receivers have shown the organization of gravity 
waves/TIDs traveling equatorward from the poles. Figure 28 shows TIDs evident in simul- 
taneous GPS and airglow data over Japan, clearly revealing the large scale density variations 
associated with such events (Saito et al. 2001). TIDs are generally grouped into large scales 
(>1000 km) and mesoscales, (scalelengths between 300 and 1000 km), the latter being 
termed MSTIDs in the literature. The fact that such waves include electric fields which can 
easily map along the magnetic field lines explain conjugate MSTID measurements observed 
by Shiokawa et al. (2004) and Otsuka et al. (2004). The coupling processes between the E 
and F region ionosphere are now being modeled (e.g., Yokoyama et al. 2009) and under- 
score the importance of gravity wave sources in understanding the mid-latitude ionospheric 
structure. 
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Fig. 29 The shear in the neutral 
wind creates a layer of enhanced 
ionization (after Hines 1964) 



Using measurements from the Dynamics Explorer-2 satellite gathered near the end of 
its mission between 250-300 km, Earle et al. (2008) reported gravity wave measurements 
at mid-latitudes, comparing simultaneous waveform measurements of the neutral density, 
plasma density, ion velocity, and neutral velocity perturbations. They showed plasma density 
variations (AN/N) of approximately 30 % with scale lengths along the satellite orbit of 
over 100 km, and argued that because the neutral density fluctuations are observed to lag 
their neutral vertical velocity by 90 degrees (phase quadrature), the perturbations are gravity 
waves as predicted by Hines (1964). They also showed that the observations agree with 
gravity wave dissipation models of Vadas and Fritts (2005) and Vadas (2007). These authors 
argued that such measurements are more common at the lower altitudes, below where most 
ionospheric satellite probes generally gather their measurements. 

4.6 Distinct Metallic Ion Layers in the Lower Ionosphere (“Sporadic-E” Layers) 

In Sect. 2, we showed that the main ions in the lower ionosphere (<150 km) were NO + 
and O 2 ■ Missing from this description of the near-earth plasma environment is the fact that 
this region is often characterized by distinct, heavy ion layers that are common in both the 
daytime and nighttime ionosphere and may be somewhat more prevalent at mid-latitudes. 
Such layers are often termed “sporadic-E” layers, since they occur in the lower E-region 
and may abruptly influence radiowave applications as they drift within the field-of-view 
of beams of HF radiowave communication antenna. In fact, the layers are actually much 
more ubiquitous than this name implies. We nevertheless retain this nomenclature for the 
discussion here. 

Sporadic-E layers are dramatic examples of neutral-plasma coupling that routinely oc- 
cur in the earth’s lower ionosphere/thermosphere system. Typically forming in the altitude 
region of 90-130 km, it is fairly well established that such layers are set up through a com- 
bination of neutral wind shear and the effects of long-lasting metallic ions resulting from 
meteoric input (e.g., Whitehead 1970; Mathews 1998; Haldoupis 2011, this issue). That 
sharp, distinct sporadic-E layers result from the occurrence of shears in the zonal neutral 
wind, as advanced by Hines (1964), Whitehead (1970), and others, is shown conceptually 
in Fig. 29 (after Hines 1964). Such layers may then be modulated by upward propagating 
gravity waves and shear instabilities (see discussion below), possibly creating polarization 
electric fields within the crests and troughs of the modulated layer. 
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Fig. 30 Incoherent scatter measurements of the plasma density in the lower ionosphere for two consecutive 
days at the Arecibo Observatory in Puerto Rico. Note the sporadic-E layer continues unbroken between day 
and night 


To illustrate the character of sporadic-E layers, altitude profiles of continuous plasma 
density observations from the Arecibo incoherent scatter radar are shown in Fig. 30 for two 
full days, in a format similar to that of Fig. 9, except the focus is now on the E-region. 
The narrow layers of enhanced ionization that are more distinct at night but are also present 
during the day are the sporadic-E layers. Of particular importance is the fact that these layers 
transition from day to night without any apparent break, and they descend, in general, with 
local time, although occasionally the layers “jump” to slightly higher altitudes. 

In addition to research carried out with radars and ionosondes, sounding rocket experi- 
ments have provided considerable insight to the formation and character of sporadic-E lay- 
ers. In particular, large amplitude neutral winds and neutral wind shears have been observed 
in vapor trails released along near vertical rocket trajectories and shown to exist at the same 
altitudes of simultaneously observed sporadic-E layers (see, for example, Larsen et al. 1998, 
and references therein). Those authors also showed calculations of the Richardson number 
based on the measured wind profile, demonstrating that the altitude range near the layer was 
highly unstable. 

In another experiment, a sounding rocket was launched in conjunction with observa- 
tions of sporadic-E layers recorded by the incoherent scatter radar at Arecibo, Puerto Rico. 
The top panel of Fig. 31 presents sporadic-E layers between 100 and 110 km in continuous 
plasma density profiles obtained with the Arecibo radar during the evening when a sounding 
rocket was launched from the north coast of Puerto Rico at 21:45 LT, as indicated by the 
arrow. The rocket payload included an ion mass spectrometer from which data in the lower 
panel of Fig. 3 1 clearly shows that the lower density layer was composed primarily of Fe + 
and Mg + ions (Earle et al. 2000). These results are in agreement with previous mass spec- 
trometer measurements of metallic ions associated with sporadic-E layers (see Grebowsky 
and Aikin 2002, and references therein). Earle et al. (2000) used the observed separation 
with altitude of the different metallic layers to estimate the zonal winds necessary to sustain 
the layer against diffusive effects. 

The combined action of neutral wind shears, electric fields, and the sharp plasma density 
gradients intrinsic to the sporadic-E layer subsequently drive instabilities that include fea- 
tures that are detectable by radars on the ground. In fact, an important consequence of the 
wind shears is the fact that they can be unstable to the Kelvin-Helmholtz instability. Such 
a process has been shown by Larsen (2000) to account for the quasi-periodic (QP) radar 
echoes that have been regularly reported in backscatter radar data associated with sporadic-E 
layers (e.g., Yamamoto et al. 1991; Hysell et al. 2004). The effects of the Kelvin-Helmholtz 
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Arecibo Observatory 
Incoherent Scatter Plasma Density 
24-25 March 1998 



Launch 


NASA 21.115 — Tortuguero Bay, Puerto Rico 
24 March 1998 at 01:45:21 U.T. 



Fig. 31 Incoherent scatter measurements of the plasma density on the night of March 24-25, 1998 at the 
Arecibo Observatory in Puerto Rico ( upper panel). A sounding rocket was launched at 01:45:21 UT (see 
arrow in top panel ) which recorded the ion composition data shown in the lower panel (after Earle et al. 
2000) 


instability on a low altitude, nearly horizontal vapor trail released from a sounding rocket 
in the presence of sporadic-E and QP-echoes have been analyzed by Larsen et al. (2005). 
In this experiment, periodic whirls are clearly visible in the vapor emissions emitted along 
a horizontal trajectory in the E-region, suggesting that the Kelvin-Helmholtz instability was 
at work during this experiment. 
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Fig. 32 Electric field measurements recorded over an active thunderstorm by two rockets, a balloon, and a 
ground receiver (after Kelley et al. 1985). The sketch on the left shows the experiment geometry. The panels 
on the right show data from each platform for one half-second period when one rocket was in the ionosphere 
and one in the mesosphere. Note how the DC electric field pulse is observed simultaneously in both rocket 
data sets and appears to penetrate upwards, into the ionosphere. The large amplitude, very low frequency 
sinusoidal variations in the upper two panels are due to the spin of the rockets 


Finally, we mention another layer phenomenon commonly observed in the mid-latitude 
ionosphere below roughly 200 km — the so-called “intermediate” or “descending” layer 
(e.g., Earle et al. 1998; Mathews 1998). An example is shown in the Arecibo density pro- 
files in the top panel of Fig. 3 1 . A descending layer appears to evolve from a density bulge 
associated with the lower edge of the F-region near 22 FT, from which a low amplitude, but 
distinct, density enhancement transitions to lower altitudes during the next hour or so. The 
creation and dynamics of such layers are not fully understood. 

4.7 Electrical Contributions to the Upper Atmosphere Due to Tropospheric Storms 

The lower edge of the ionosphere has been recognized for some time as the conductive 
upper layer of the earth-ionosphere waveguide. As such, it provides an essential part of 
the global circuit driven by atmospheric electricity (e.g., Wilson 1920). The ability of low 
frequency (DC-EEF) electric fields associated with tropospheric storms to penetrate into the 
ionospheric plasma remains a topic of considerable discussion (see Simoes et al. 2011, this 
issue). 

To investigate the extent of electrical energy deposition in the ionosphere, Kelley et al. 
(1985) reported observations gathered by a pair of sounding rockets launched simultane- 
ously to apogees in the mesosphere and ionosphere over an active thunderstorm at a mid- 
latitude location (Wallops Island, Virginia, USA). As shown in Fig. 32, the results clearly 
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demonstrated that the DC electric field “pulse” associated with lightning spherics penetrated 
into the ionosphere. Besides requiring a revision of models of the transfer of electromag- 
netic energy across the atmosphere/plasma interface, these results demonstrated that sig- 
nificant energy is deposited in the ionosphere by lightning-generated electric fields. These 
researchers also reported a significant component of the DC electric field pulse was oriented 
parallel to the magnetic field direction in the lower ionosphere. Such an electric field would 
be expected to contribute to localized plasma heating, as discussed by those authors. 

The effects of tropospheric storms on the ionosphere may include other consequences. 
Burke et al. (1992) reported Dynamics Explorer-2 observations of Alfven waves near 
350 km altitude over a hurricane that show a DC electric field pulse of more than 40 mV/m 
with an associated plasma density depletion and strong magnetic signature. These authors 
argued that this signature was due to Alfven “wings” associated with an intense upwards- 
propagating DC electric field generated by the hurricane forces. Similarly, using data from 
the San Marco satellite, Farrell et al. (1994) reported large, localized DC electric field “struc- 
tures” in the ionosphere over tropical storms, providing further evidence that DC-coupled 
electric fields associated with tropospheric storms propagate upwards and penetrate into the 
ionosphere. These and numerous other observations underscore how the ionosphere is sig- 
nificantly influenced by energy and momentum sources originating in the troposphere below. 
We now turn to the high latitude ionosphere and consider energy and momentum sources 
that originate from above, in the magnetosphere. 


5 The High Latitude Ionosphere 

The ionosphere at high latitudes differs significantly from that at low and mid lati- 
tudes. At high latitudes, the ionosphere exchanges energy and momentum with the solar 
wind/magnetosphere primarily via the highly inclined magnetic field lines that characterize 
the earth’s polar regions. The magnetic field thus serves as a conduit for both charged parti- 
cles and their associated field-aligned currents and DC electric fields, thus directly coupling 
the high latitude ionosphere and upper atmosphere with the magnetosphere. Although photo- 
ionization still contributes to the ambient plasma environment, an alternative, often domi- 
nant, source of ionization in the high latitude ionosphere is impact ionization of energetic 
particles with the neutral atmosphere, and the primary drivers of high latitude ionosphere 
dynamics are electric fields that map down from the magnetosphere. In contrast, the pri- 
mary source of ionization for the low and mid-latitude ionosphere is photo-ionization from 
the sun and the primary drivers of plasma transport are neutral winds and small amplitude 
electric fields created by the neutral winds. 

In this section, we introduce the visual aurora and the different populations of precipi- 
tating electrons that create the aurora. We discuss the creation of thermal plasma by auroral 
precipitation and compare this with photoionization. We then review the dynamics and cur- 
rents of the high latitude ionosphere associated with magnetospheric processes. This is fol- 
lowed by a discussion of the many fundamental processes in which the high latitude plasma 
environment interacts with the upper atmosphere. We conclude with some brief remarks 
regarding auroral arc electrodynamics. 

5.1 The Visual Aurora and Energetic, Precipitating Auroral Particles 

Overview of the Visual Aurora. Optical emissions associated with the aurora result from in- 
teractions of energetic, precipitating particles with the upper neutral atmosphere. Typically 
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Fig. 33 Photograph of the aurora taken in January, 2005 at Bear Lake, Alaska. Credit: Joshua Strang, USAF 



Fig. 34 View of the aurora photographed from the International Space Station on March 4, 2012, 
17:20:33 UT on a pass over the South Indian Ocean. Figure courtesy NASA 


produced within the altitude region of 80-250 km, displays of discrete auroral arcs are per- 
haps the most striking, visible manifestations of the interactions of space plasma — in this 
case, precipitating energetic electrons — with atmospheric, neutral gases. Figure 33 is an au- 
roral display photographed from the ground that shows sheets of auroral emissions, which 
are modulated by large amplitude spatial undulations that have typical durations of seconds 
to minutes. Another view of such auroral “curtains” is shown in Fig. 34. In this case, the 
side of the aurora has been photographed from the NASA International Space Station, at an 
altitude near 400 km. Note that the auroral emissions are characterized by nearly vertical 
striations that correspond to the ambient magnetic held lines along which the precipitating 
particles are guided. The green glow at the base of the aurora is the most common auroral 
hue (A. = 557.7 nm) and is due to interactions of the more energetic electrons with oxygen 
at the lower altitudes (100-130 km) whereas the pink-red glow (A = 630.0 nm) in the upper 
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portion of the auroral display is from interactions of less energetic electrons with oxygen 
at higher altitudes (130-250 km). Other, less common, colors of auroral displays include 
bright red and blue, which correspond to emission lines associated with the interaction of 
precipitating energetic particles with other gases such as N 2 . 

An image of the aurora gathered by the DMSP satellite above the southern polar region 
is shown in Fig. 35 and depicts discrete arcs, diffuse aurora, and auroral forms filled with 
turbulent-like structure. This image demonstrates that, in the aggregate, auroral emissions 
and the processes that create them are essentially global in their extent, expanding across 
the polar region with a variety of patterns in which shorter scale structures are embedded. 
The ensemble often resembles an auroral oval (sometimes referred to as an auroral belt), as 
is illustrated in the images with the wider field of view shown in Fig. 36. This time sequence 
of images, gathered every 12 minutes by the imager on the Dynamics Explorer- 1 satellite, 
illustrates how dynamic the large scale aurora may be, in this case brightening and expanding 
poleward on the nightside during a magnetospheric substorm. (See Frank and Craven 1988 
for a discussion of the auroral images from this instrument.) 

The source of the auroral optical emissions shown in the preceding figures is energetic 
electrons with typical energies of 100’s of eV to 10’s of keV that precipitate into the high 
latitude neutral gases where they impact and are “braked”. Here, such energetic electrons 
interact with the neutral atmosphere in two basic ways: excitation and ionization. Excitation 
occurs when an energetic electron strikes a neutral atom or molecule and excites the neu- 
tral to a higher energy state. When the neutral relaxes to a lower energy state, a photon is 
released, creating the visible aurora. Ionization of the ambient neutral gas occurs when the 
incoming energetic electrons impact the neutral atoms or molecules and liberate electrons 
to create ions, thermal electrons, and more energetic secondary electrons. In the remainder 
of this sub-section, we discuss the different types of aurora and their relative distributions in 
the high latitude ionosphere. We then discuss the ionization associated with the precipitating 
particles in the following sub-section. 

The aurora may be categorized into three primary types that depend on their respective 
acceleration or precipitation processes and on the characteristics of their associated pre- 
cipitating energetic particles: (1) discrete, (2) oscillating or Alfvenic, and (3) diffuse. Each 
are introduced below, followed by comparisons of their relative number flux and energy 
flux contributions to the high latitude ionosphere. Although we mention other auroral fea- 
tures as they pertain to the near-earth plasma environment, we do not attempt to review all 
of the aurora-related processes here. Readers are directed to many reviews of auroral phe- 
nomena and related physical processes, including those by Rees (1989), McFadden et al. 
(1999), Lummerzheim (2007), and Newell et al. (2009). In particular, the ISSI monograph 
entitled, Auroral Plasma Physics, edited by Paschmann et al. (2003), contains a thorough 
compendium of current understanding of the physical processes associated with the aurora. 

Discrete Aurora. The visual aurora shown in Figs. 33, 34, 35, 36 are, in general, as- 
sociated with discrete auroral arcs. In situ measurements by probes on sounding rock- 
ets and satellites have shown that the precipitating energetic electrons associated with 
discrete auroral arcs have narrow, nearly monoenergetic spectra (e.g., Mcllwain 1960; 
Evans 1968). These electron spectra are typically organized in “inverted- V”-shaped energy 
functions that correspond to spatial variations across the discrete auroral arc (Frank and 
Ackerson 1971). (The inverted-V’s tend to resemble inverted-U’s where the energy is dis- 
played as a logarithm.) The peak energies of such energetic electrons typically range from 
100’s of eV to a few tens of keV. In Fig. 37, energetic electrons detected near 4000 km 
altitude by electrostatic analyzers on the FAST satellite are shown alongside an image of 
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Fig. 36 Consecutive false-color images of the northern auroral oval recorded at the onset and evolution 
of a substorm as the Dynamics Explorer 1 spacecraft descends from near its apogee during the period 0529 
through 0755 UT on 2 April 1982. Each image is telemetered in a twelve minute interval. The poleward bulge 
of an auroral “surge” is first seen in the fourth image, which begins at 0605 UT. A westward-traveling surge 
then propagates along the poleward edge of the oval as highly structured, eastward-moving auroral forms 
develop in the post-midnight sector. The sunlit atmosphere is visible at the upper left in each image. The 
passband of the optical filter is 123-155 nm, in which the dominant responses are from emissions of atomic 
oxygen at about 130.4 and 135.6 nm and the LBH bands of molecular nitrogen (figure courtesy NASA. See 
also Frank and Craven 1988) 


multiple auroral arcs simultaneously photographed by an all-sky camera (Stenbaek-Nielsen 
et al. 1998). Although the optical emissions are from much lower altitudes (near 100 km), 
the correspondence of the visual arcs with the in situ measurements of the higher altitude 
electron structures is evident. 

Significant progress has been made in the last 15 years establishing the relation be- 
tween accelerated, energetic electrons and the discrete aurora, primarily using observations 
gathered from satellites that included probes with high time resolution designed specifi- 
cally to address these processes, such as Viking, Freja, and FAST. As discussed by nu- 
merous authors (e.g., Carlson et al. 1998; McFadden et al. 1999; Paschmann et al. 2003; 
and Marklund 2009), data from these satellites have helped establish that electrostatic “U- 
shaped” potentials are set up in what is often referred to as the auroral acceleration re- 
gion, existing at altitudes of roughly 2,000-20,000 km, high above where the electrons 
interact with the upper atmosphere and the auroral light is emitted. Such potential struc- 
tures contain electric field components oriented along the magnetic field direction that 
account for the nearly monoenergetic electrons responsible for the discrete aurora. The 
same potential structures also accelerate ion beams along the magnetic field in the op- 
posite direction, as observed by these same satellite probes (e.g., Shelley et al. 1976; 
McFadden et al. 1999). 

An important feature of the data in Fig. 37 is that on the equatorward side of the event 
near 9:18:15 UT, there is an inverted-V region of upwards accelerated electrons. These elec- 
trons have no atmosphere in which to interact and thus no visual analog in the all-sky 
camera display. Indeed, the lack of optical emissions in such circumstances has been re- 
ferred to as the “black aurora” by Marklund et al. (1994). High resolution satellite probes 
on the Viking, Freja, and FAST satellites (e.g., Marklund et al. 1994; Carlson et al. 1998; 
McFadden et al. 1999) have identified such upwards accelerated electrons as a regular 
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Fig. 37 Multiple arcs in a 
low-light “all sky” image 
recorded on an aircraft {top) and 
by instruments on the FAST 
satellite {lower panels) near 
21:20 MLT. In the top panel, the 
FAST trajectory is mapped along 
the magnetic held line to the 
100 km conjugate point at 10 s 
intervals as the satellite passes 
across from left to right. The 
second panel is the precipitated 
energy flux on a linear scale. The 
narrow auroral arcs are clearly 
displayed here, as are their 
correspondence to the visible 
aurora in the top panel. The 
auroral optical emissions, in 
particular the two arcs to the far 
right (north), changed over the 
4 minutes of the pass. The two 
color panels show the measured 
electron energy spectra for 
downgoing (third panel from 
bottom) and upgoing (second 
panel from the bottom) and show 
a number of inverted-V 
structures. Note the strong 
upgoing electrons present in the 
FAST data at the start of the 
event near 09 : 1 8 : 1 5 UT that have 
no visual analog in the all-sky 
image. The lower panel shows 
the FAST zonal magnetometer 
data corresponding to the 
downgoing current region, where 
the upgoing electrons are 
observed, and the upgoing 
current region, where the 
downgoing electrons are 
observed (after Stenbaek-Nielsen 
et al. 1998) 
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feature of the aurora. The satellite probes have also measured the perpendicular electric 
fields in the V- shapes and have shown how they reverse polarity in accordance with the 
oppositely-directed accelerated particles (e.g., Marklund et al. 1994; Carlson et al. 1998; 
McFadden et al. 1999). Although such upwards-directed electrons are accelerated at high 
altitudes, they are believed associated with current closure in the lower ionosphere, where 
the current is carried by lower energy electrons. The manner in which the lower ionosphere 
controls this current closure remains one of the fundamental outstanding problems in auroral 
physics today. 

Oscillating Auroral Electrons Associated with Alfven Waves. Another means by which 
electrons are accelerated by electric fields and subsequently precipitate into the upper atmo- 
sphere to create aurora is their interactions with Alfven waves. In this case, the oscillating 
electric fields in the wave accelerate field-aligned electrons in both directions along the 
magnetic field which are often referred to as “counter-streaming” accelerated electrons. An 
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Fig. 38 Observations from the FAST satellite on the polar cap boundary. Top: the electric field along the 
spacecraft trajectory, or roughly north-south. Middle: the energetic electron spectra. Bottom: the pitch angle 
spectra of the energetic electrons with 0° pointing downwards (after Chaston et al. 2004) 


example gathered by the FAST satellite is shown in Fig. 38 (Chaston et al. 2004). Notice that 
the pitch angle information shows the electrons as highly aligned along, and anti-parallel to, 
the magnetic field direction, although the energies themselves are spread-out over a large 
range of energy space, typically between 10 eV and 1 keV. Other detailed observations of 
Alfven-wave driven field-aligned electrons, often within an inverted-V event characterized 
by higher energy electrons, have been reported in sounding rocket data sets (e.g., Chen et al. 
2005). 

Although not originally appreciated as an important source of aurora, recent observations 
by probes on the FAST satellite have shown that such Alfven wave accelerated electrons can 
account for as much as 1/3 of the total energy deposited in the ionosphere during elevated 
auroral activity (Chaston et al. 2007). Wygant et al. (2000) showed how the Poynting flux 
associated with strong electric and magnetic signatures of Alfven waves observed on the 
high-altitude Polar satellite corresponded to, and was sufficient to power, intense auroral 
optical emissions. 

Diffuse and Pulsating Aurora. Diffuse aurorae are considerably fainter than the vi- 
sual aurora associated with discrete auroral arcs and are best discerned with the aid 
of photometers. Although such aurorae are generally without distinct, sharp features, 
there is some structure that can be discerned from the measurements (Lui et al. 1973; 
Johnstone et al. 1993). In general, the diffuse aurora is most prevalent in the post-midnight 
sector and has wider latitudinal extent than the discrete aurora. 

The diffuse aurora may include large scale pulsations in which large patches of aurora 
faintly turn on and off with periods typically lasting from 1-30 seconds. Correlated with 
the recovery phase of magnetic sub-storms, the pulsating aurora is also prevalent in the 
post-midnight hours and, as with the diffuse aurora, is believed caused by wave-particle 
interactions with whistler mode waves (e.g., Johnstone 1983; Davidson 1990; Nishimura 
et al. 2010). Recent observations suggest that ionosphere-magnetosphere coupling may play 
an important role in the generation mechanism (Sato et al. 2004). 
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DMSP 11-year Energetic Electron (32eV - 30 keV) Distributions - Number Flux 
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Fig. 39 Comparisons of average number fluxes of the discrete aurora (left), Alfvenic aurora (middle), and 
diffuse electron aurora (right), for low solar wind driving conditions (top) and high solar wind driving con- 
ditions (bottom). The data were recorded by the energetic particle instrument on the DMSP satellites over 
the interval of 1988 to 1998. The total number flux for each group is shown (after Newell et al. 2009). The 
outermost circle in each case corresponds to 50° and corrects an error in the original manuscript (Newell, 
personal communication, 2012) 


Unlike the discrete and Alfvenic aurorae which are produced by energetic electrons ac- 
celerated by components of electric fields along the magnetic field direction, diffuse aurora 
are caused by pitch angle scattering with whistler mode chorus waves in the magnetosphere 
(Thorne et al. 2010), and then precipitate when their loss cones intersect the upper atmo- 
sphere. Despite the fact that the diffuse aurora does not create as strong an optical signature, 
it deposits considerably more energy than the other forms of aurora, as we discuss next. 
Diffuse ion aurora also exists, and precipitates via the same mechanism, as discussed, for 
example, by Newell et al. (2009). 

Relative Energy Flux of Different Auroral Processes. Newell et al. (2009) have averaged 
energetic electron data from the polar-orbiting DMSP satellites over an 1 1 year time from 
1988 to 1998. The data were binned according to the type of auroral acceleration (discrete, 
Alfvenic, or diffuse) and the results are shown in Figs. 39, 40 corresponding to averages of 
number flux (Fig. 39) and energy flux (Fig. 40). The data in each case were further binned 
by low and high solar wind driving conditions. As explained by those authors, the local time 
coverage of the satellites was fairly uniform, although the post-midnight region was covered 
only from 1992-1998. 

The distributions shown in Figs. 39, 40 compiled by Newell et al. (2009) immediately 
show that discrete aurora are most prevalent on the dayside and pre-midnight portions of 
the auroral oval, depositing most of their energy in the post-noon, pre-midnight sector. The 
Alfvenic aurora is a dayside and midnight phenomenon, yet deposits most of its energy on 
the nightside, just prior to midnight. The diffuse aurora is primarily a post-midnight/morning 
phenomenon, although it is also the most prevalent form of aurora in the cusp centered at 
noon on the dayside. Of paramount importance is that the diffuse electron aurora is clearly 
the dominant source of particle energy flux into the high latitude ionosphere, constitut- 
ing 63 % of the energy in the low solar wind driving conditions and 57 % during the 
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DMSP 1 1 -year Energetic Electron (32eV - 30 keV) Distributions -- Energy Flux 
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Fig. 40 Comparisons of average energy fluxes of the discrete aurora {left), Alfvenic aurora (middle), and 
diffuse electron aurora (right) for low solar wind driving conditions (top) and high solar wind driving con- 
ditions (bottom). The data were recorded by the energetic particle instrument on the DMSP satellites over 
the interval of 1988 to 1998. The total energy flux for each group is shown (after Newell et al. 2009). The 
outermost circle in each case corresponds to 50° and corrects an error in the original manuscript (Newell, 
personal communication, 2012) 


high solar wind driving conditions, according to Newell et al. (2009). Note that the ion 
diffuse aurora were also included in the statistics presented by these authors but are not 
shown here. That the diffuse electron aurora is so prevalent has important consequences for 
the conductivity of the high latitude, post midnight ionosphere, as we discuss further on 
below. 

Cusp Aurora. The cusps are funnel-shaped magnetic field regions in the northern and 
southern hemispheres on the dayside of the earth directly connected with the solar wind. 
As such, they permit direct entry of solar wind particles and energy into the high latitude 
ionosphere/upper atmosphere. The cusp aurora thus represents a distinct region of the near- 
earth plasma at high latitudes. Indeed, the ion diffuse aurora peaks in the cusp (Newell et al. 
2009). 

An example of data gathered by probes on the FAST satellite as it crossed the cusp 
region near its perigee at 348 km is shown in Fig. 41. The lower panel shows the down- 
going, energetic ions of solar wind origin, exhibiting an energy dispersion indicative of a 
cusp crossing in the presence of the southward interplanetary magnetic field (e.g., Reiff 
et al. 1977). The precipitating electrons are organized in very narrow, inverted-V forms 
with modest peak energies of 100’s of eV and thus deposit their energy and plasma 
at higher altitudes, in the 200-300 km region (as discussed in the next sub-section). 
The presence of such inverted-V signatures suggest that local electron acceleration pro- 
cesses in the cusp are similar to those in the auroral acceleration region on the nightside 
(e.g., Yamauchi et al. 1998; Pfaff et al. 1998). The upper panel shows the highly corre- 
lated electric and magnetic field signatures, indicative of field-aligned currents and down- 
ward directed Poynting Flux present in the cusp regions (e.g., Strangeway et al. 2000; 
Le et al. 2002), as we will discuss later. 
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Fig. 41 Observations by instruments on the FAST satellite recorded in the cusp region near 348 km altitude 
and local magnetic noon. The bottom and middle panels show energy spectra of highly structured downgoing 
ions and electrons, respectively. The top panel shows the electric field component along the spacecraft velocity 
(red) and the zonal component of the magnetic field (black), which are highly correlated 


Fig. 42 Thermal plasma 
production rates (left) and 
resulting density profiles (right) 
as a function of incoming 
electron energy (Strickland et al. 
1983) 
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5.2 Thermal Plasma Created by Precipitating Energetic Electrons 

The interaction altitude of the energetic electrons with the atmosphere depends on both the 
energy of the precipitating electrons and on the neutral density and composition of the atmo- 
sphere. The effective penetration altitude of precipitating electrons as a function of energy 
is shown in Fig. 42 for a model atmosphere (Strickland et al. 1983). Alongside is a panel 
that shows the ion production rates associated with these interactions. As with the mid and 
low-latitude plasma discussed previously, the molecular ions (NO + and 07) produced at the 
lower altitudes (100-130 km) recombine quickly. At higher altitudes (i.e., 200-300 km), the 
plasma created by the softer electrons (i.e., a few hundred eV) that ionize the 0 + , have much 
longer lifetimes and hence are important in subsequent transport and dynamic coupling to 
the neutral gas, as we discuss later on below. 

Another example of precipitating, energetic electrons measured by the FAST satellite 
at 437 km altitude and organized in inverted-V shaped energy functions with peak energies 
ranging from 100’s of eV to several keV is shown in the upper panel of Fig. 43. As discussed 
above, the inverted-V shapes represent spatial variations of the electron energy that occur 
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FAST 5 December 1998 
ALT = 437 km MLT=15.6 hours 
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Fig. 43 Observations recorded on the FAST satellite near 437 km are shown in the top panel and reveal 
narrow beams of downgoing, energetic electron spectra. The bottom panel shows the corresponding thermal 
plasma created by these electrons, according to the electron transport model described in Strickland et al. 
(1993) 

across auroral arcs. The lower panel shows the expected thermal plasma associated with 
these measured electron energies, according to the electron transport model described in 
Strickland et al. (1993). Note how the altitude of the model thermal plasma in the lower 
panel varies in accordance with the measured energy of the precipitating electrons. 

Incoherent scatter radar observations reveal how thermal plasma created by keV electrons 
impinging on the atmosphere is indeed extended in altitude within the lower ionosphere. 
Figure 44 (figure courtesy I. Haggstrom) shows a time sequence of continuous data gathered 
with a fixed radar beam directed up the magnetic field line, and hence the same space-time 
cautions that were discussed with respect to the radar data presented in Fig. 24 must again be 
invoked here. These radar observations were gathered simultaneously with all-sky camera 
observations showing a series of “intense, high, rayed arcs”, as reported by Kirkwood et al. 
(1988). 

Diffuse aurora, which is characterized by precipitation of electrons with a broad energy 
spectrum (from 100’s of eV to tens of keV) over a wide latitudinal extent, can last for 
hours, creating substantial ionization over a wide range of altitudes (e.g., 100-300 km). The 
ionization at the higher altitudes persists much longer due the much longer recombination 
times of the 0 + ions, as noted above and discussed in Sect. 2. 

5.3 Characterizing the High Latitude Ionosphere Plasma 

The high latitude plasma density is directly formed via: (1) precipitating energetic electrons 
that create thermal plasma as they are braked by the atmosphere (as discussed above) and 
(2) photo-ionization. We next discuss the contribution from photo-ionization at high lati- 
tudes. In subsequent sections, we discuss how other processes must be taken into account, 
particularly transport via E x B drifts and winds, in order to specify the state of the high 
latitude ionosphere at any given location and time. 

Photo-Ionization at High Latitudes. In the preceding sub-sections we showed how the 
aurora is statistically distributed within the high latitude region and discussed how the re- 
sulting thermal plasma varies as a function of altitude, depending primarily on the energy of 
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Fig. 44 EISCAT incoherent scatter radar observations of plasma density recorded with a beam oriented in a 
fixed position up, along the magnetic field direction (figure courtesy I. Haggstrom) 


the precipitating electrons and the background neutral density and composition. Although 
ionization produced by such energetic auroral particles often dominates the local high lati- 
tude ionosphere, photo-ionization at high latitudes is nevertheless an extremely important, 
regular source of plasma in these regions as well. 

Contributions to the plasma density from photo-ionization in the high latitude ionosphere 
calculated using the IRI model for summer and winter hemispheres at the December solstice 
are shown in Fig. 45. Although the solar elevation angle is low at high latitudes, the contin- 
uous sunlight during the summer creates a conductive E-region over the polar cap to 70 de- 
grees on the midnight side (lower panel). During the winter, on the other hand, the opposite 
conditions occur, and the E-region is devoid of plasma from the midnight sector over the 
poles to 70 degrees on the noon side (upper panel). As with the mid- and low-latitude iono- 
sphere, recombination efficiently diminishes the plasma population below about 200 km in 
the absence of sunlight, yet the F-region ionosphere is maintained, albeit at lower densities, 
even during the winter months at the pole. 

Plasma Temperature and the Polar Wind. The plasma temperature in the high latitude 
ionosphere is highly variable due to the many different and highly dynamic processes at 
work in this region. Precipitating energetic electrons heat the ambient plasma to higher tem- 
peratures, although the temperature of the secondary electrons as well as heating contri- 
butions from associated local plasma waves remain controversial, as discussed in Wahlund 
et al. (1993). Those authors used EISCAT radar data to show that electron temperature en- 
hancements associated with 100-500 eV electron precipitation correlated well with plasma 
density enhancements in the altitude region of 170-230 km. Heating via other mechanisms 


Springer 





The Near-Earth Plasma Environment 


69 


Polar Cap Electron Density — IRI Model (2007) 
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Fig. 45 IRI model calculations (for December 22, 2004) of electron density in the ionosphere starting at 
50 degrees latitude, crossing the pole, and extending to 50 degrees latitude on the other side. Top: north- 
ern hemisphere. Bottom: southern. Evidence for numerical deviations in the model at the geographic pole 
precludes accurate estimates in this region 



such as frictional Joule heating driven by perpendicular electric fields (discussed below), 
and plasma instabilities also create regions of increased plasma temperatures (e.g., Jones 
et al. 1992; Williams et al. 1992). Although this heating is initially localized, the increased 
pressure gradients subsequently affect neighboring regions. Despite the strong variability of 
the auroral precipitation within the high latitude ionosphere, the electron temperature has 
been modeled for different altitudes of this region by Schunk et al. (1986). 
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A consequence of the larger temperatures and the nearly vertical magnetic held direction 
in the polar region is that plasma hows outward along the magnetic held lines, often with 
signihcant velocities. This “polar wind” is a feature of the polar ionosphere that results 
from an ambipolar electric held set up as the electrons rise to higher altitudes along the 
magnetic held lines due to their higher temperature and lighter mass, as discussed in reviews 
by Ganguli (1996) and Schunk and Sojka (1997). The polar wind includes both 0 + and H 1 
ions that eventually escape the ionosphere into the magnetosphere and distant tail regions. 
Observations of the polar wind have been measured by instruments on numerous spacecraft, 
including Dynamics Explorer-1 (Chappell et al. 1982, 1987; Nagai et al. 1984), Akebono 
(Abe et al. 1993), and Polar (Su et al. 1998). It is important to note that other forms of ion 
outhow, in which ions are accelerated to higher energies via wave-particle interactions, also 
help populate the magnetosphere with ions of ionospheric origin, as reviewed by Andre and 
Yau (1997) and others. 

Structure in the High Latitude Ionosphere. The high latitude ionosphere is highly struc- 
tured, which results from a variety of mechanisms too numerous to review here. (See, Kelley 
2009, for a wide-ranging discussion.) The most obvious source of plasma structure is pre- 
cipitating auroral particles, as illustrated in the auroral images in Fig. 33, 34, 35 and 36. 
These structures can involve scales sizes from 1000’s of km to less than a km. Such density 
forms have well-defined plasma gradients and associated electric field and current structures 
that drive local instabilities and produce irregularities with associated scintillations. 

An important feature of the structured high latitude ionosphere is the appearance of 
patches of enhanced plasma density in the F-region that convect over the polar cap in the 
anti-sunward direction. These patches have horizontal scales that range from 100-1000 km 
and have associated smaller scale structure and irregularities that also produce radio scin- 
tillation effects (e.g., Tsunoda 1988; Coley and Heelis 1995). The patches appear when 
the interplanetary magnetic field (IMF) is southward and have been shown by numerous 
researchers (e.g., Buchau et al. 1985; Valladares et al. 1996) to result not from energetic 
particle precipitation, but rather from detached dayside F-region plasma that is subsequently 
transported over the polar cap as part of the general anti-sunward flow pattern (discussed 
below). This narrow, sunward stream of enhanced dayside aurora is referred to as a “tongue 
of ionization” and is strongest when magnetic activity is high (e.g., Sojka et al. 1981; 
Middleton et al. 2005; Foster et al. 2005). 

The ionospheric plasma is particularly structured within the cusp region. Using a survey 
of 18 months of EFF electric field data from the Dynamics Explorer-2 satellite, Heppner 
et al. (1993) showed that the most intense low frequency (4-8 Hz) structures at high latitudes 
clearly are those in the dayside cusp region. These authors interpreted their electric held 
variations as either static structures or Alfven waves. Measurements gathered along vertical 
profiles provided by a sounding rocket, as shown in Fig. 46, show that the plasma density is 
highly structured along the portion of the trajectory within the cusp proper on open magnetic 
held lines, yet is smooth on that portion of the trajectory that is outside of the cusp measured 
on closed magnetic held lines. These data were gathered in conjunction with low frequency 
wave electric held irregularities and are interpreted primarily as static structures that created 
strong scintillations that were observed on the ground during the flight (Pfaff et al. 2003). 

5.4 Farge-Scale Electrodynamics of the High Fatitude Ionosphere: Magnetospheric 
Electric Fields and Field- Aligned Currents 

We now discuss the electrodynamics of the high latitude ionosphere, for which we focus on 
the coupling of this plasma to the magnetosphere via DC electric helds and field-aligned 


Springer 


The Near-Earth Plasma Environment 


71 


Fig. 46 Plasma density along 
vertical profiles, measured by a 
sounding rocket in the cusp 
( lower panels). The upper panel 
shows how the rocket payload 
was launched in the cusp region 
from Ny-Alesund and remained 
there during the upleg portion of 
the trajectory. The rocket then 
crossed from open to closed 
magnetic fields lines during the 
downleg portion of the trajectory, 
where the plasma irregularities 
and structure are no longer 
observed 



Downleg Upleg 



currents. Later, we introduce effects of the neutral atmosphere. Space limitations restrict our 
discussion here to descriptions of general features relevant to the dynamical and electrical 
coupling of the ionosphere and magnetosphere. For more detailed discussions, readers are 
referred to numerous publications on this subject, including Lyons (1992), Richmond and 
Thayer (2000), Cowley (2000), and Paschmann et al. (2003). 

Overview of the Electrical Connection Between the Solar Wind, Magnetosphere, and 
Ionosphere. The solar wind is a magnetized plasma that flows around the earth’s mag- 
netic envelope, or magnetosphere. Figure 47 depicts the magnetosphere in the noon- 
midnight/north-south plane in which the earth is shown to scale in the center. For reference, 
the radius of the earth is ~6400 km. Although the near-earth plasma region that is the subject 
of this article would only extend a short distance on this diagram, the earth’s magnetic field 
enables this region to communicate and exchange energy and momentum with the farthest 
reaches of the magnetosphere. 

Since the solar wind is supersonic, a bow shock and underlying magnetosheath region 
is formed upstream of the compressed magnetosphere which then deflects impinging solar 
wind particles. The outer boundary of the magnetosphere is a current layer called the mag- 
netopause and exists where the solar wind pressure and the earth’s magnetic field pressure 
are equal. On the sunward side, the magnetopause is typically located at 8-12 earth radii 
(Re), depending on the solar wind dynamic pressure. 
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Fig. 47 The main characteristics of the magnetosphere and its interaction with a B- < 0 interplanetary 
magnetic field. The dawn-dusk electric field is shown, along with the principal direction of the plasma flow. 
The toroidal volume near the earth ( shaded orange) is the co-rotating plasmasphere (adapted from Lyons 
1992) 

Solar wind energy connects through the magnetopause to the magnetosphere via its mag- 
netic field. The presentation shown in Fig. 47 depicts the standard Dungey-type “X-line” 
configuration of interaction of the earth’s magnetosphere with the solar wind with the 
IMF B- < 0 (Dungey 1961). (The standard coordinates used for describing the sun- 
magnetosphere system are: x, directed from the center of the earth towards the sun, y, 
directed from dawn to dusk, and z, directed upwards, perpendicular to the ecliptic.) The 
X-line represents a merging of the magnetic fields of the earth and of the solar wind and 
is referred to as magnetic “reconnection” (e.g., Priest and Forbes 2000). In cases where the 
IMF B- < 0, significant energy is transferred from the solar wind to the magnetosphere, 
whereas, when IMF B z > 0, this transfer is reduced. 

An essential feature of solar wind/magnetosphere interaction described above is that a 
dawn-dusk electric field is set up as a result of the solar wind flow across the earth’s magnetic 
field dipole, given by E = — V.vw x B. This electric field creates anti-sunward E x B plasma 
flow in the polar cap regions. The polar cap regions are essentially characterized by “open” 
magnetic field lines that do not close in the conjugate hemisphere of the earth, but are swept 
back to form the earth’s magnetotail. Equatorward of these open field lines is the plasma 
sheet, which is a region of trapped plasma on closed magnetic field lines that are stretched 
down the tail and yet whose “feet” map to the auroral oval regions of the high latitude 
ionosphere in each hemisphere. This region corresponds to the annular regions where auroral 
optical emissions are most common, as shown earlier in the images in Fig. 36. The direction 
of the electric field within this ring in the dawn-dusk direction is opposite that in the polar 
cap, which is key to understanding plasma flows in the high latitude ionosphere. 

To illustrate this important feature, consider the magnetosphere depicted from the per- 
spective of looking from the magnetotail towards the sun in Fig. 48 in which only the upper 
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Fig. 48 The magnetosphere from the perspective of the tail, looking toward the sun. The principal current 
systems and electric fields are shown. The Pedersen currents with which the magnetospheric currents close 
in the ionosphere are generally much closer to the earth than shown (adapted from Lyons 1992) 


half of the magnetosphere is shown. In the polar cap region above the earth at the center of 
the illustration where the magnetic field lines are essentially vertical, the dawn-dusk electric 
field extends directly to the earth’s lower ionosphere. On the flanks of the polar cap adjacent 
to the magnetopause, the return plasma flow is sunward at the equator and the associated 
duskward-directed electric field on each side maps back (earthward) along the magnetic 
field lines to create an electric field oriented in the opposite direction to that in the polar 
cap. In this manner, we see that the electric field associated with the return, sunward flow, 
when mapped to regions closer to the earth, is oriented towards the pole on the dusk side 
and towards the equator on the dawn side. In each case, these electric fields drive sunward 
plasma flow in the auroral oval. 

The electric field orientation and resulting plasma flow in the high latitude ionosphere 
is now shown from a polar perspective, in the left hand side of Fig. 49(a). The polar cap 
and auroral oval electric fields combine to drive a pair of vortices, or convection cells, 
that are responsible for the plasma circulation within the entire high latitude ionosphere. 
The shapes and orientations of the dual cell patterns are highly dependent on forcing by 
the interplanetary magnetic field (IMF), as discussed by many authors (e.g., Heelis 1988; 
Cowley 2000). 

In addition to the patterns of high latitude electric fields, field-aligned currents of magne- 
tospheric origin are also organized in a highly systematic manner. The principal field-aligned 
currents that connect the magnetosphere and ionosphere are those that flow at the boundaries 
of the auroral oval in each hemisphere. These currents are depicted in Fig. 49(b) where their 
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High Latitude Ionosphere 
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Fig. 49 The electrodynamics of the polar cap region for an idealized, symmetric coupling to the magneto- 
sphere. The illustration on the left (a) shows the electric fields and E x B drifts; the illustration on the right 
(b) shows the principal current systems (inspired by Prolss 2004) 


relation to the two-celled convection pattern and high latitude electric fields is apparent. The 
field-aligned currents on the poleward side (the so-called “region 1” currents) connect to the 
magnetopause, as shown in Fig. 48. Driven by stresses on the dayside magnetopause (e.g., 
Siscoe and Siebert 2006), the region 1 currents flow out of the ionosphere on the duskside 
and into the ionosphere on the dawn side. These currents then flow as Pedersen currents in 
the ionosphere and return to the magnetosphere on the equatorward side of the auroral zone, 
where they are referred to as “region 2” currents, flowing into the ionosphere on the dusk 
side and out of the ionosphere on the dawn side. 

A second set of currents is driven within the high latitude ionosphere by the applied 
electric fields discussed above. These fields drive not only Flail and Pedersen currents, given 
by (6) in Sect. 3 and also shown in Fig. 49(b), but may also drive field aligned currents in 
the following way. If the conductivity of the ionosphere were uniform, the E x B plasma 
flows might exist without the need for a field aligned current system to maintain equilibrium. 
This is not the case in nature. Gradients in the conductivities force ionospheric currents to 
flow along field lines. Furthermore, as pointed out by Cowley (2000), the Pedersen currents 
flowing along the electric field direction in the ionosphere cannot close in the ionosphere, 
but rather a system of field-aligned currents is needed to accommodate their divergence. 
Note that these field-aligned currents are in addition to the principal, region- 1 field-aligned 
currents driven by the stresses at the dayside magnetopause. 

We now present experimental evidence that supports the main features of the electrical 
coupling between the high latitude ionosphere and magnetosphere described above. All of 
the data shown were gathered by probes flown on polar-orbiting satellites in the high latitude 
ionosphere. 

High Latitude Convection. Measurements of DC electric fields and plasma drifts have 
been fundamental in establishing the plasma circulation in the high-latitude ionosphere. 
Among the earliest satellite measurements of the high latitude DC electric field are those 
from the OGO-6 satellite (Heppner 1972), for which an example is shown in Fig. 50. These 
data were gathered on a dusk-to-dawn pass of the satellite across the polar ionosphere. The 
large amplitude electric fields in the auroral oval which drive sunward flow are clearly vis- 
ible. The polar cap electric fields are considerably weaker, but clearly correspond to anti- 
sunward E x B flow. 
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Fig. 50 An example of the horizontal DC electric held component perpendicular to the sun-earth line, ob- 
served along the dusk-dawn meridian from a low altitude pass of the OGO-6 satellite over the northern polar 
cap (Heppner 1972) 


Direct measurements of the DC electric fields and plasma drifts in the high latitude 
ionosphere by numerous authors including Heelis (1984), Heppner and Maynard (1987), 
and others have shown how high latitude convection changes in strength and pattern as 
a function of the IMF. In the case where the IMF B z < 0, the plasma drifts are not only 
stronger, but they also show a much more organized structure with two well-defined cells 
that may be skewed dawnward or duskward according to the direction of the IMF By com- 
ponent. In cases where IMF B~ is >0, the convection is not well organized, and the ob- 
served fields (drifts) are weaker. In fact, during such weakly driven conditions, the potential 
patterns are such that additional, isolated convection cells may be present (Heelis et al. 
1986). 

Models of the high latitude convection patterns, such as those based on ionospheric elec- 
tric fields measured by the Dynamics Explorer-2 satellite (Weimer 1995) show how the 
two-celled patterns respond to the interplanetary magnetic field, as shown in Fig. 51. This 
particular presentation shows the potential patterns, which correspond to plasma convection 
flow lines. These patterns underscore how the general flow of the high latitude ionosphere 
is governed by the interplanetary magnetic field. 

Field-Aligned Currents. In concert with large-scale patterns of magnetospheric electric 
fields, measurements of field-aligned currents have helped to confirm the basic magneto- 
spheric energy flow processes discussed above. The charges carried by the precipitating 
energetic particles that create the aurora discussed in previous sub-sections have an associ- 
ated current that flows along the direction of the magnetic field. Although such field-aligned 
currents were first proposed by Birkeland (1908) with respect to the aurora, their existence 
was not fully accepted until space-based magnetometer measurements showed their pres- 
ence definitively (Zmuda 1966; Cummings and Dessler 1967). 

The large-scale field-aligned currents that flow between the ionosphere and magneto- 
sphere are not necessarily associated with optical aurora displays. In fact, these currents, 
which are almost always present, may be carried by low energy electrons that would not 
necessarily have sufficient energy to create the impact ionization or optical emissions dis- 
cussed above. By convention, the field-aligned current is designated upward (downward) for 
downgoing (upgoing) electrons. 

In situ magnetometer data from polar-orbiting satellites revealed that such field-aligned 
currents are typically observed as a pair of infinite sheets, providing important clues re- 
garding how the majority of the currents might close. The upper portion of Fig. 52 depicts 
field-aligned currents as flowing vertically along an infinite current sheet that is traversed 
in the normal direction by a spacecraft carrying a vector magnetometer. The left hand il- 
lustration of a single, upward current would produce a step function in the east-west (or 
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Fig. 51 Models of the high-latitude electric potential as a function of the interplanetary magnetic field (IMF), 
using spherical harmonic coefficients derived from measurements of the electric held by the Dynamics Ex- 
plorer-2 satellite (Weimer 1995). The models shown here are for cases in which the IMF magnitude, Bj , was 
greater than 7.25 nT. See Weimer (1995) for an explanation of details on the figure panels 




zonal) component of the satellite magnetometer data. The right hand side shows the mag- 
netic perturbation from two anti-parallel sheets. Magnetometer data from the MAGSAT 
satellite gathered along a dusk-to-dawn polar orbit are shown in the lower portion of the 
figure (Zanetti et al. 1984). The observed zonal magnetometer component is consistent 
with two parallel current sheets of oppositely directed vertical currents encountered on each 
side of the polar cap. Note that the observed magnetometer signature appears of opposite 
sign when converted to east-west coordinates on the equatorward-bound satellite trajec- 
tory. 
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Fig. 52 (Top) How data from the 
east-west component of a vector 
magnetometer along a polar orbit 
can be used to deduce 
field-aligned currents in single 
and dual parallel current sheets 
(Pfaff 1996). (Bottom) Fluxgate 
magnetometer data recorded on a 
polar crossing by the MAGSAT 
satellite. The measurements show 
two regions of parallel current 
sheets, with field-aligned currents 
into and away from the 
ionosphere, encountered on each 
side of the polar cap (Zanetti 
et al. 1984) 
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Statistical studies of these events showed that the field-aligned currents form a character- 
istic pattern, arranged in a belt around the polar cap corresponding to where auroral optical 
emissions are customarily observed, as shown in Fig. 53 (Iijima and Potemra 1976). Further- 
more, the patterns show that the currents on the poleward edge of the polar cap flow into the 
ionosphere on the dawn side and out on the dusk side, with opposite flows on the equator- 
ward edge of the oval. These measurements were fundamental to demonstrating the patterns 
of Region 1 /Region 2 currents about the auroral oval that connect the ionosphere and mag- 
netosphere, as shown in Fig. 48. Similar patterns have been recorded in other satellite data, 
for example in the FAST data set by Peria et al. (2000). 

In cases where the field-aligned currents are carried primarily by accelerated energetic 
particles, detailed comparisons of data such as those gathered on the FAST satellite show 
that zonal magnetometer perturbations match extremely well the field- aligned currents com- 
puted from the measured fluxes of field-aligned, energetic electrons (Elphic et al. 1998). For 
example, returning to Fig. 41 that shows FAST satellite observations within the cusp, notice 
that in each instance where downgoing packets of electrons were observed, the sign of the 
zonal magnetic field component is negative, indicative of upgoing currents. Such measur- 
able variations, even for extremely small electron fluxes as shown in this figure, illustrate 
how commonplace field-aligned currents are in the earth’s high latitude plasma. 

The primary manner in which field-aligned currents are believed to “close” in the iono- 
sphere is by horizontal Pedersen currents, flowing in the direction of the ambient DC electric 
field, as sketched in Fig. 48 and Fig. 49(b). Although it is expected that the majority of such 
currents would close via a neighboring current sheet as discussed above, these illustrations 
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Fig. 53 A summary of the 
distribution and flow directions of 
large-scale, field-aligned currents 
determined from data obtained 
on 493 passes of the Triad 
satellite during weakly disturbed 
conditions (|AL| < lOOy). The 
hatched area shown between 
1130 and 1230 MLT in the polar 
cusp region indicates the current 
flow directions are not distinct 
(Iijima and Potemra 1976) 
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t- "i Currents Away from Ionosphere 


show that some of the field-aligned currents might close via Pedersen currents over the polar 
cap. In fact, Le et al. (2010), using ST-5 satellite data, demonstrated that such polar cap Ped- 
ersen currents were balanced by an observed asymmetry in the measured Region 1/Region 2 
currents. 

Closure of the field-aligned currents occurs in the lower ionosphere, where the Pedersen 
conductivity is largest. It is controlled both by the variations with altitude in the local ther- 
mal plasma (i.e., local conductivity), determined primarily by the spectra of precipitating 
energetic electrons, and by the Pedersen mobilities that depend on the charged-neutral col- 
lision frequencies and thus the properties of the atmospheric gas (see Sect. 3.2). The sunlit 
ionosphere on the dayside provides increased conductivity via photo-ionization to facilitate 
this closure. 

Correlation of Electric and Magnetic Fields and Poynting Flux. Observations of DC elec- 
tric fields in the high latitude ionosphere show a very strong correlation with the magnetic 
field variations, as shown in Fig. 54 (Sugiura et al. 1982). The strong correlation is generally 
present between the poleward component of the electric field, which we designate as A E x , 
and the eastward component of the perturbation magnetic field (i.e., in which the earth’s 
ambient magnetic field has been subtracted), which we designate as A B y . These observa- 
tions support the infinite current sheet approximation for the field- aligned currents, and is 
consistent with meridional closure, i.e., via Pedersen currents. Using this assumption and 
assuming that the Pedersen conductivity is uniform along the measurement track, the two 
quantities may be related by AB y /AE x = /x c S p , where F p is the height integrated Pedersen 
conductivity and /x 0 is the permeability of free space (e.g., Smiddy et al. 1980). 

Such combined electric field and magnetic field measurements are also important for cal- 
culating the Poynting flux into the ionosphere (e.g., Kelley et al. 1991). The Poynting flux is 
given by AE x AB and its divergence can be related directly to the energy input to the iono- 
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Fig. 54 Data from the east-west 
magnetic field component (after 
subtraction of the earth’s ambient 
magnetic field) and the 
north-south electric field 
component measured on the 
polar-orbiting Dynamics 
Explorer-2 satellite (Sugiura 
et al. 1982) 
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Fig. 55 Ion drift velocity ( upper panel ) and magnetic field ( middle panel) components in the anti-sunward 
direction recorded by the Dynamics Explorer-2 satellite for a pass over the northern polar region. The bottom 
panel shows the field-aligned Poynting flux S|| = S ■ B c , which is primarily from these two components. The 
data are from Gary et al. (1994), which shows all measured components 


sphere that is dissipated as Joule heat, as discussed further in the next section. In Fig. 55, 
we present an example of field-aligned Poynting flux, S||, measured by instruments on the 
Dynamics Explorer-2 satellite (Gary et al. 1994). In this presentation, the ion drift compo- 
nent is shown, which can be equated to the E x B component at these altitudes (>350 km). 
The data shown here were gathered on a dusk-to-dawn traversal of the polar cap in which 
the components of the ion drift and magnetic field were measured in the anti-sunward direc- 


<£) Springer 


80 


R.F. Pfaff 


tion. (See Gary et al. 1994 for a detailed discussion as well as a presentation of all 6 of the 
measured components.) Not only does this hgure again show the correlation of the crossed 
electric and magnetic held components, but it emphasizes how the electromagnetic energy 
is directly delivered to the high latitude ionosphere and upper atmosphere — i.e., downward, 
along the magnetic held direction. 

Alfven Waves. Alfven waves are essentially ubiquitous in the high latitude ionosphere, as 
the medium is not in a steady state and motions of “slabs” of plasma (e.g., enhanced conduc- 
tivity regions corresponding to auroral arcs), or shears in the electric held, will launch Alfven 
waves which will carry field-aligned currents (e.g., Kelley 2009). Furthermore, Alfven “res- 
onators” have been proposed as regions between large gradients in the plasma density, such 
as between the F-region and the lower portion of the magnetosphere near (say) 3000 km 
(e.g., Lysak 1993), where Alfven waves may be trapped. Note that such resonators need not 
be conhned to the high latitude region. 

From an experimental standpoint, following the discussion above concerning the high 
degree of correlation between electric and magnetic helds, we note that the electric and 
magnetic held signatures within Alfven waves are also highly correlated and can mimic 
“static structures”, particularly since the static structures can be easily Doppler-shifted by 
the satellite motion into the same frequency regime as the Alfven waves. An examination of 
the ratios of the measured A E and A B components can be compared to the expected Alfven 
velocity, Va, which is given by V a = |5 g |/(4jt / o) 1 / 2 , where, p is the particle density which 
depends on the plasma number density and the ambient ion composition. A comparison of 
estimates of the Alfven velocity and the integrated Pedersen conductivity (discussed above) 
can then be used to distinguish the origin of the correlated electric and magnetic helds (e.g., 
Knudsen et al. 1990). Such comparisons using in situ measurements have been discussed in 
detail by Gurnett et al. (1984), Berthelier et al. (1989), Knudsen et al. (1990), Heppner et al. 
(1993), and others. 

5.5 Atmospheric Interactions with the Magnetospheric Energy and Momentum 

In this sub-section, we introduce the role of the neutral atmosphere within the framework of 
the interactions of the high latitude ionosphere/upper atmosphere with the magnetosphere. 
We have already shown how the atmosphere brakes incoming charged particles, creating 
thermal plasma as well as auroral optical emissions, and have noted the atmosphere’s role 
with respect to current closure of field-aligned currents. We now address additional pro- 
cesses that involve the high latitude upper atmosphere and its interactions with magneto- 
spheric electric fields and particles: 

1. E x B Drift/Atmosphere Interactions 

- Ion drift deceleration below 250 km 

- Winds driven by magnetospheric electric helds 

- Joule heating 

- Electric helds driven by winds and upwards Poynting hux 

- Two-stream waves driven by electrons streaming through collision-dominated ions 

2. Energetic Particle/Atmosphere Interactions 

- Creation of thermal plasma by energetic particle precipitation 

- Heating due to particle precipitation 

- Formation of Hall and Pedersen currents (in conjunction with the electric held) and 
held- aligned current closure 
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Fig. 56 Ion drift velocity 
perpendicular to the magnetic 
field in the zonal {red) and 
meridional {black) directions, 
and the E x B drift in the zonal 
direction {blue), as a function of 
altitude for a DC electric field of 
50 mV/m applied in the 
meridional directional and 
assumed to be constant with 
altitude. The variations of the ion 
drift component with altitude are 
largely due to ion-neutral 
collisions, as explained in the 
text. The slight variation of the 
E x B magnitude with altitude is 
due to small changes of the 
magnetic field strength with 
altitude 


Poker Flat, Alaska 


E-50 mV/m meridional 



Ion Drift Deceleration Below 250 km. As discussed above, the DC electric fields that 
map down from the magnetosphere set the high latitude plasma in motion via E x B drifts. 
At lower altitudes (below 250 km), however, the ion motion is appreciably affected by ion- 
neutral collisions such that the ion plasma drift is noticeably slowed and undergoes a change 
in direction. The electrons, on the other hand, remain fully “magnetized” and continue to ex- 
ecute E x B drifts at altitudes as low as 80-90 km. Below these altitudes, the electron-neutral 
collision frequency becomes significant when compared to the electron gyro frequency, as 
shown earlier in Sect. 3. 

The ion and electron plasma drifts, calculated versus altitude for a 50 mV/m poleward 
(meridional) DC electric field at Poker Flat, Alaska (|B| = 0.54 G), are plotted in Fig. 56. 
Ions and electrons both undergo an E x B drift of approximately 950 m/s at 250 km. How- 
ever, at lower altitudes, where the ion-neutral collisions become increasingly substantial, the 
ion drift begins to both slow and change its orientation towards the direction of the DC elec- 
tric field. At about 125 km altitude, where the Pedersen mobility peaks, the ion drift is one 
half the value given by E x B and is equally divided between flow along and perpendicular 
to the electric field direction. As the ion-neutral collisions slow the ion drifts, momentum is 
transferred to the neutral gas, as discussed next. 

Winds Driven by the Magnetospheric Convection. The strong plasma motions in the high 
latitude ionosphere that are driven by the magnetospheric DC electric fields set the neutral 
upper atmosphere in motion via ion-neutral collisions. Thus, the electromagnetic energy of 
the magnetosphere is converted to mechanical energy. To illustrate this coupling, Fig. 57 
shows simultaneous measurements of the ion drifts (associated with the E x B drift) and the 
neutral wind measured by the Dynamics Explorer-2 satellite as it transited the polar region 
near 350 km altitude (Killeen et al. 1984). The neutral winds essentially track the ion drifts in 
the evening sector and within the polar cap, yet diverge considerably in the morning sector. 
Where the two vectors are not aligned is where Joule heating and compositional changes 
occur, as discussed below. 
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Fig. 57 The neutral wind and 
ion drift vectors measured by 
Dynamics Explorer-2 on a 
perigee pass over the southern 
polar region, plotted in 
geographic polar coordinates. 
The neutral winds are coded by 
the yellow arrows, and the ion 
drifts are coded by the red bars. 
The cur\'ed line represents the 
location of the solar terminator 
(90° solar zenith angle). The 
symbol S refers to the location of 
the invariant pole at the universal 
time of the given pass (Killeen 
et al. 1984) 



As the high latitude neutral motions are influenced by both the amplitude and direction 
of the E x B drifts associated with the magnetospheric electric fields, the resulting neu- 
tral winds show a strong and remarkable IMF dependence in accordance with the E x B 
convection, as discussed by Ki lleen and Roble (1988) (see, in particular, their Fig. 18). 

Although the neutral wind circulation is expected to follow that of the plasma convection 
at altitudes of 300 km and above, since the ion velocities decrease and rotate with decreas- 
ing altitude below 250 km (see Fig. 56), the neutral wind circulation at lower altitudes is not 
well established. The situation is further complicated by the fact that the neutral winds are 
also driven by tidal forcing from below. Furthermore, during auroral events, the plasma con- 
ductance is often increased within certain altitude intervals corresponding to thermal plasma 
created by energetic particles, altering the coupling efficiency of the ions and neutrals. For 
cases where steady plasma channels might be created via diffuse auroral precipitation, atmo- 
spheric models predict that neutral “jets” may form (Brinkman et al. 1995). The existence of 
such jets and the neutral circulation at low altitudes associated with magnetospheric electric 
fields are largely unexplored in any systematic way. 

Joule Heating. The relative motions of the ion and neutral drifts have a profound effect 
on the temperature of the neutral gas via Joule heating. Joule heating may be considered 
the result of friction between the ion and neutral gases that occurs when the gases are not 
flowing at the same speed and direction. The Joule heating rate per unit volume, Qj, may 
be represented as: 

Qj = °> E ' 2 (?) 

where E' = (E + UxB) and a p is the Pedersen conductivity, discussed in Sect. 3.2. Clearly, 
the Joule heating depends linearly on the plasma density and on the square of the relative 
motions of the plasma drifts and neutral winds. Joule heating has a strong altitude depen- 
dence due to the variation of the Pedersen conductivity and of the winds with altitude. The 
neutral gas heating is also affected by the fact that cooling rates increase with decreasing 
altitude as the neutral density increases. 

To illustrate local Joule heating, we return to the example in Fig. 57, but now present in 
Fig. 58 other quantitative measurements gathered in situ by instruments on the Dynamics 
Explorer-2 satellite during this same event. The neutral winds track the ion drifts throughout 
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Fig. 58 Geophysical quantities 
measured along the track of 
Dynamics Explorer-2 during the 
orbit shown in Fig. 20(a). The ion 
drifts and the neutral winds are 
shown in the top two traces, 
plotted against time, altitude, and 
latitude of the spacecraft. The 
second panel shows the electron, 
ion, and neutral temperatures 
measured along the track, and the 
third panel shows the atomic 
oxygen and molecular nitrogen 
number densities ( left-hand 
scale) and the electron density 
(right-hand scale). The bottom 
trace shows the ion-neutral 
coupling time constant measured 
along the track, as discussed in 
Killeen et al. (1984) 
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most of the satellite pass over the polar cap prior to 23:36 UT and at these times, there is little 
local Joule heating expected and hence there is little change in observed ion temperature. 
Notice, however, how the ion temperatures increase at later times when the drifts diverge 
(Killeen et al. 1984). This is also a region where compositional changes are observed caused 
by thermal expansion. 

It is important to note that structured electric fields will produce significantly more Joule 
heating than smoothed fields that typically are used in model calculations (Codrescu and 
Fuller-Rowell 1995). Indeed, various researchers (e.g., Thayer 1998; Emery et al. 1999) have 
invoked short scale electric fields to account for a “missing energy” in circulation models 
and other related calculations. Joule heating represents a fundamental process whereby mag- 
netospheric electromagnetic energy is deposited into the upper atmosphere in the form of 
heat, and, assuming a steady state and motionless neutral gas, can be shown to be equivalent 
to the Poynting flux (e.g., Thayer and Semeter 2004). Joule dissipation, frictional heating, 
and the related energy exchange in the ionosphere/thermosphere system have been discussed 
by numerous authors (e.g., Thayer et al. 1995; Thayer 1998; 2000; Strangeway 2012) and 
represent an important means by which the atmosphere acts as a reactive load on the mag- 
netosphere. 

Wind Driven Electric Fields and Upwards Poynting Flux at High Fatitudes. The winds 
set up by the plasma drifts driven by magnetospheric electric fields have other profound con- 
sequences in the high latitude regions. The winds, unbound by the magnetic field, transport 
energy and momentum outside the regions of the plasma drifts. Furthermore, due to their 
inertia, they do not readily change direction as the applied E x B direction changes nor do 
they cease to flow when the driving DC electric fields are turned off. In other words, the high 
latitude neutral winds have persistent scales and time constants that are longer than those 
corresponding to the applied electric fields. Consider the case where the neutral gases are 
set in motion and then the driving magnetospheric electric fields are shut off abruptly. In this 
case, the neutral gases maintain their proper motion and generate an associated electric field 
that then acts back on the system. This “flywheel” effect then generates Poynting flux di- 
rected upwards, into the magnetosphere, opposite to the customary Poynting flux associated 
with the impinging magnetospheric energy (Kelley et al. 1991). Such “upwards” Poynting 
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Fig. 59 Altitudes in the high 
latitude, lower ionosphere 
(Kiruna, Sweden) where the 
two-stream instability is unstable 
according to the linear dispersion 
relation, as a function of the 
driving DC electric field (Pfaff 
1986) 



flux has been observed on the Dynamics Explorer-2 spacecraft by Gary et al. (1994) during 
IMF B z > 0 conditions. 

Two-Stream Waves in the High Latitude, Lower Ionosphere Regions. The collisional two- 
stream instability represents another common mechanism in the high-latitude lower iono- 
sphere in which energy is extracted from the moving plasma driven by magnetospheric 
electric fields. This process occurs in the lower ionosphere, where E x B driven electrons 
“stream” through the ions that are collisionally bound to the neutral gas. A two-stream insta- 
bility occurs when the relative ion-electron velocity exceeds a threshold given approximately 
by the ion acoustic velocity, as discussed earlier with respect to the equatorial electrojet. In 
the auroral electrojet, such relative ion-electron drifts occur at altitudes between 95 and 
1 25 km where the Hall current is dominant. For typical magnetic field strengths at these al- 
titudes within the high latitude regions, the instability threshold is met for DC electric fields 
of approximately 20-25 mV/m, since the ion acoustic velocity here is about 350-400 m/sec. 
Although the two-stream instability is typically associated with the auroral electrojets, this 
condition may be met throughout the high latitude ionosphere, including the polar cap, pro- 
vided the DC electric field is of sufficient magnitude. 

The altitudes where the two-stream instability is generally met are shown in Fig. 59 as 
calculated for Kiruna, Sweden, by Pfaff (1986). The unstable region is given by the thresh- 
old condition of fluid theory (e.g., Sudan et al. 1973) given by V, — V e = C s ( 1 + \//). Here, 
C s is the acoustic velocity and is the ratio of the product of the ion-neutral and electron- 
neutral collision frequencies to the product of the ion and electron gyro frequencies. Al- 
though ifr is about 0.1 at 105 km at high latitudes and decreases for higher altitudes, this 
variable grows significantly larger at lower altitudes and accounts for the sharp cut off of the 
instability at altitudes of about 95 km and below, as shown in the figure. At the higher alti- 
tudes, the instability shuts off as the ions become magnetized and the relative ion-electron 
velocity is reduced. Electric field observations of two stream waves measured on a sound- 
ing rocket launched from Esrange, Sweden are shown in Fig. 60 and are described in Pfaff 
et al. (1992). In general, the waves correspond to meter scales and typically display large 
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Fig. 60 Observations of DC 
electric field magnitude (left) and 
wave electric fields driven by the 
two-stream instability (right) 
measured on the downleg of a 
sounding rocket flight from 
Kiruna, Sweden. These data are 
described in Pfaff et al. (1992) 


Flight 21 .097 -- Esrange, Sweden 
9 April 1988 
Downleg 



DC Electric Fields 
125- 

120 - 

115- 

£ HO' 

-S 105- 

s 100- 

95- 
90- 

85 J T r- 

0 50 100 

mV/m 


125 

120 

115 

110 

105 

100 

95 

90 

85 


AC Electric Fields 


1000 2000 3000 4000 
Frequency (Hz) 


-10 

r>j 



O 

O 


-50 


amplitudes in which the wave electric fields are a significant fraction of the driving DC 
electric field. Furthermore, the waves are known to heat the local electrons, raising the elec- 
tron temperature by factors of 2-4 (e.g., Schlegel and St.-Maurice 1981; Jones et al. 1991; 
Williams et al. 1992). Although driven by DC electric fields set up by an entirely different 
process, the same instability is also common in the equatorial electrojet, as discussed in 
Sect. 4. 

Creation of Thermal Plasma and Heating by Energetic Particle Precipitation. In 
Sect. 5.2, we discussed impact ionization and how the interactions of energetic particles 
of sufficient energy create thermal plasma as they are braked by the upper atmosphere. This 
mechanism represents a fundamental means by which the atmospheric gas and the magne- 
tospheric plasma interact. 

Absorbing the kinetic energy of incoming particles heats the neutral gas. On average, 
as reported by Lu et al. (1995), this heat source is about | that of typical Joule heating, 
although it can be significant during large auroral events (e.g., Thayer and Semeter 2004). 
The heating caused by the mechanical energy of the precipitating particles is important for 
the high-latitude energy budget as well for the creation of localized effects. For example, 
Clemmons et al. (2008) have invoked heating by soft particle precipitation in the cusp to 
suggest that upwellings might explain the depletions in the neutral density measured by the 
Streak satellite at low altitudes. 

Formation of Hall and Pedersen Currents and Field-Aligned Current Closure. In the 
lower ionosphere, where the ion drifts are slowed by collisions, differential drift motion 
of the ions and electrons drives a Hall current. As discussed with respect to the equato- 
rial electrojet, this Hall current maximizes at the lower altitudes where the ions become 
essentially collision dominated. As shown in Fig. 49(b), on the nightside, the pre-midnight 
poleward DC electric field drives the electrons westward in the auroral oval, which results 
in an eastward electrojet, and the opposite is true for the equatorward DC electric field at 
post midnight, where a westward electrojet is driven. The electrojet currents close in the 
polar cap via the weaker Hall currents, which depend on both the available plasma density 
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as well as the driving electric fields, and thus the current path is not straightforward, par- 
ticularly where both of these quantities are highly variable. Where the currents converge in 
the midnight sector, known as the Harang Discontinuity, as well as in the cusp region on 
the dayside, field-aligned currents may be produced to keep the overall currents divergence 
free. 

Pedersen currents are also driven by the applied electric fields, as discussed above with 
respect to the closure of the main Region 1/Region 2 field-aligned currents. In discussing the 
Pedersen currents it is important to realize that whereas the Pedersen mobility peaks near 
125 km, the altitude profile of the thermal plasma on the nightside auroral zone is highly 
dependent on the energy of the incoming energetic particles, as shown earlier in Fig. 43, and 
hence the Pedersen conductivity may peak at a different altitude than the Pedersen mobility. 
In other words, as the Pedersen conductivity is dependent on the combined mobility and 
plasma density effects, we see that the efficiency of the closure currents for a given auro- 
ral precipitation event is thus intricately tied to both the energy spectrum of its energetic 
particles and the characteristics of the ambient atmosphere. 

5.6 Auroral Arc Electrodynamics and Feedback to the Magnetosphere 

We conclude this section by briefly returning to the subject of discrete auroral arcs and their 
associated electrodynamics. Given the relevance of the aurora to the physics of the high 
latitude ionosphere, a few comments regarding their inherent physical processes as well as 
their role in coupling the ionosphere and magnetosphere are appropriate. 

We showed previously how energetic precipitating auroral particles impact the neutral 
atmosphere to form narrow channels of enhanced thermal plasma within a limited altitude 
range, generally below 300 km. In many ways, the creation of thermal plasma is simply 
one step in a chain of processes that includes the acceleration of the energetic particles 
and the feedback of the enhanced conductivity on the generation mechanism itself. The 
many related electrodynamic processes associated with discrete auroral arcs include: (1) the 
closure of the current associated with the precipitating particles, (2) the effect of the newly 
formed auroral arc plasma on the ambient DC electric field, (3) the enhancement of electric 
fields within the arc and along its sides, and (4) the feedback of the enhanced conductivity 
on the acceleration mechanism responsible for accelerating energetic particles. The reader 
is referred to the various chapters in the ISSI monograph edited by Paschmann et al. (2003) 
for a detailed discussion of many of these predominantly plasma processes. 

Among the earliest observational results from in situ measurements gathered in the 
auroral ionosphere by sounding rockets is that of the anti-correlation of DC electric 
fields and plasma density within auroral arcs (e.g., Aggson 1969; Potter and Cahill 1969; 
Evans et al. 1977; Maynard et al. 1973; Mallinckrodt and Carlson 1985). To first order, 
such an anti-correlation is explained via Ohm’s law, J = a ■ E. In order to maintain current 
continuity, the electric held is reduced where the conductivity (i.e., inside the arc) is in- 
creased. How the modified electric held effects the source potential in the magnetosphere 
remains an open question, particularly in cases where the conductivity of the entire au- 
roral oval is signihcantly enhanced (e.g., Fig. 36). The picture is further complicated by 
the motions of the enhanced conducting “slabs” which affect electric helds on neighbor- 
ing magnetic held lines. Several researchers have shown that auroral arcs have a proper 
motion that is different from that of the ambient E x B drifts (e.g., Haerendel et al. 1993; 
Frey et al. 1996). Indeed, Mallinckrodt and Carlson (1985) suggested that a larger collisional 
coupling with the neutrals may be necessary to explain the different velocities, and obtained 
a relationship between the neutral wind and the cross-arc electric held potential. 
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With respect to the electrodynamics within the arc, the observed electric fields are not 
always decreased (e.g., Marklund 1984; Paschmann et al. 2003). The creation of polar- 
ization electric fields within auroral arcs has been explained by invoking enhanced Cowl- 
ing conductivities within the arc channel (e.g., Bostrom 1973; Paschmann et al. 2003). 
Other researchers have reported enhanced electric fields and associated velocity shears at 
the edges of auroral arcs (e.g., Kelley and Carlson 1977). Measurements with the EIS- 
CAT radar have shown that the plasma drifts at the edges of the arcs may be in op- 
posite directions (e.g., Haerendel et al. 1996), which may eventually help explain au- 
roral motions such as curls and folds (e.g., Hallinan and Davis 1970; Hallinan 1976; 
Davis 1978). 

As discussed earlier, it is now fairly well established that electric field structures (some- 
times referred to as “shocks”) at altitudes of several thousand km and higher accelerate 
charged particles along magnetic field lines, and that electrons accelerated downward by 
these structures impact the atmosphere and create the discrete aurora. These electric field 
structures are set up to enable the ionospheric plasma to maintain continuity of the field- 
aligned currents. If the available conductivity in the volume in which the currents must flow 
is too small, as explained by Knight (1973), electric fields oriented parallel to the magnetic 
field are generated to accelerate the plasma along the magnetic field direction and hence 
greatly increase the available current in a given volume. (See, for example, Paschmann et al. 
2003; Ergun et al. 2004; and references therein.) 

Of particular importance is the discovery of “inverse” electrostatic structures of opposite 
polarity that account for the acceleration of upgoing electrons (e.g., Marklund et al. 1994; 
Carlson et al. 1998; McFadden et al. 1999). Besides direct measurements of these acceler- 
ated particles, as shown by these authors, the satellite probes have measured the perpendicu- 
lar, U-shaped potentials that have the necessary reversed polarity to produce the oppositely- 
directed accelerated particles. The currents of these paired shocks appear to be associated 
with the field-aligned currents discussed above which close in the lower ionosphere. Further- 
more, electrodynamic models which include low altitude acceleration of upward-directed 
electrons have been advanced by Doe et al. (1995) and by Marklund et al. (1997) to explain 
what appear to be auroral density cavities at altitudes as low as 150 km in the high latitude 
ionosphere (Doe et al. 1993). These and related observations underscore our limited knowl- 
edge of how the ionosphere participates in the acceleration of energetic particles and the 
closure of field-aligned currents. 

To illustrate the importance of the ionosphere in understanding this fundamental mech- 
anism of nature, consider the role of the ionospheric conductivity in the acceleration of 
energetic particles. Newell et al. (1996) have carried out statistical analysis of aurorae using 
direct measurements of energetic particle data to demonstrate that discrete aurorae are more 
likely to be generated in darkness than in sunlight, as shown in Fig. 6 1 . These data support 
the concept that the lower conductivity of the dark ionosphere requires an auroral acceler- 
ation mechanism to be set up in order to produce the keV electrons to carry the current as 
required by the Knight relation. In a similar fashion, the abundance of the diffuse aurora in 
the post-midnight auroral zone may provide sufficient conductivity in the ionosphere such 
that the field-aligned currents are maintained without the need for potential structures to be 
set up to accelerate the particles. The diffuse aurora thus supplies the conductivity in the 
same manner as that provided by photo-ionization during sunlit conditions, in accordance 
with the fact that discrete arcs are not commonly observed during these post-midnight peri- 
ods. 

To summarize, since the ionosphere must support the field-aligned current circuit be- 
tween the ionosphere and magnetosphere, its properties are essential to maintaining these 
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SUNLIT DARK 



Fig. 61 Probability of observing discrete aurora when the polar ionosphere is sunlit {left) and dark (right) 
as a function of invariant latitude and magnetic local time, based on measurements of energetic particles 
recorded by the DMSP satellites (Newell et al. 1996) 


currents and to the formation of the electric field structures that ultimately accelerate both 
the downgoing and upgoing electrons. In this fashion, the near-earth plasma plays a critical 
role in establishing the conditions to accelerate energetic particles which create the discrete 
aurora. Because the neutral atmosphere not only brakes the incoming energetic particles but 
also enables the high conductivity thermal plasma to form and the currents to close (via 
Pedersen mobilities), the critical role that the neutral atmosphere plays in the coupling of 
the high latitude ionosphere with the magnetosphere cannot be overstated. 


6 Ionosphere Response to Magnetic Storms 

We now provide a brief discussion of some of the ways in which the near-earth plasma 
responds to solar-induced magnetic activity and in particular to magnetic storms. As we shall 
show, the entire ionosphere including middle and low latitudes undergoes dramatic changes 
during the strongest magnetic storms. As we cannot describe here all of the ways in which 
the ionosphere responds to magnetic storms, we direct the reader to a number of reviews 
on this subject, including Prolss (1995), Forbes et al. (2000), and Mendillo (2006), as well 
as the articles in the AGU monograph on mid-latitude ionospheric disturbances edited by 
Kintner et al. (2008). 

The near-earth plasma environment is continuously responding to two types of varia- 
tions of the sun: (1) variations of the flux of solar EUV radiation, as discussed in Sect. 2, and 
(2) variations of the solar wind at the earth’s location which subsequently influence the mag- 
netosphere/ionosphere/thermosphere system. In this section, we focus on the second type, 
and in particular, the response of the near-earth plasma environment to magnetic storms. 
Magnetic storms are initiated by prolonged periods of the solar wind’s southward (B z < 0) 
magnetic field direction and are most severe when associated with a “cloud” of solar wind 
plasma (e.g., Gonzalez et al. 1994). Major storms may last 2-3 days and are typically moni- 
tored by the storm-time disturbance index, or Dst. Dst is calculated from magnetometer data 
gathered at low latitude stations that provide a measure of the magnetospheric ring current 
energy. Energy transfer from the solar wind to the magnetosphere results in an increase in 
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the intensity of the ring current that corresponds to a negative Dst perturbation (e.g, Daglis 
et al. 1999). The main phase of a magnetic storm has an associated negative Dst excursion 
which lasts for several hours, with the largest storms displaying Dst values of < — 100 nT. 

The magnetic substorm is another important category of magnetic activity to which the 
near-earth plasma responds (Sergeev et al. 2012). Magnetic substorms are different than 
magnetic storms in that they are generally much briefer, and result from energy released 
from the tail of the magnetosphere which then propagates into the high latitude auroral oval. 
Substorms are characterized by brightenings and movements of auroral displays and have 
strong associated high latitude magnetic disturbances that result from increased currents. 
The appearance of substorms as well as magnetic storms are evident in an index of magnetic 
activity is known as K p that monitors worldwide magnetic disturbances every 3 -hours with 
an integer index from 0-9 corresponding to a logarithmic scale. A general discussion of 
geomagnetic indices is provided by McPherron (1995). 

The influence of magnetic storms on the near-earth plasma may be considered from two 
general perspectives: (1) magnetic storm effects that originate at high latitudes, for exam- 
ple increased magnetospheric energy and momentum deposition in the high latitude iono- 
sphere/upper atmosphere, which then have subsequent consequences at lower latitudes, and 
(2) global magnetic storm effects which are applied to the ionosphere/thermosphere system 
at all latitudes within a brief time interval after storm onset, for example those effects driven 
via prompt penetration electric fields. Below, we discuss examples of magnetic storm effects 
on the ionosphere/thermosphere system, corresponding to each category. 

6.1 Magnetic Storm Effects that Originate at High Latitudes 

Magnetospheric Energy and Particle Increases at High Latitudes. During a magnetic storm, 
the polar cap and auroral oval widen and the electromagnetic input and auroral emissions in- 
crease considerably. Although these changes are driven by the larger-scale magnetospheric 
activity, the increased energy and particle input to the ionosphere/thermosphere system is 
initially concentrated within the high latitude regions. For example, an estimate of increased 
heating at high latitudes during the magnetic storm of January 10-11, 1997, compared to 
quiet conditions, is shown in Fig. 62. This heating has been calculated by Richmond and Lu 
(2000) using a model that assimilates a variety of ground-based and space-based measure- 
ments. The upper panels show the energy input due to particle heating, whereas the lower 
panels show the energy input due to Joule heating. It is clear that the energy input to the high 
latitude region is significantly increased in both cases compared to quiet conditions and that 
this energy is greater for the Joule heating. As noted by these authors, the Joule heating has 
a greater effect on the thermospheric circulation not only because of its larger energy but 
also because the particle energy is deposited within an altitude region that extends as low as 
90 km where the neutral gas, with its significantly larger density, is more difficult to perturb. 
The storm-driven energetic particles nevertheless provide extremely important conductivity 
changes via impact ionization, as discussed in the last section, as well as other effects. 

Magnetospheric energy input at high latitudes during storms, such as shown in this ex- 
ample, can lead to significant effects at lower latitudes, as reviewed by Fuller-Rowell et al. 
(2008). For example, strong pressure gradients caused by the enhanced Joule heating in the 
high latitude regions during active periods drive strong equatorward neutral winds that are 
not constrained by the magnetic field. Neutral wind changes driven by a sudden increase 
in Joule heating at high latitudes have been computed via the CTIPe numerical simulation 
and are shown in Fig. 63 (Fuller-Rowell et al. 2008). The winds are shown at 15 UT for the 
time corresponding to three hours after the onset of the high latitude forcing. Note the wind 
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Fig. 62 Distributions in magnetic latitude and local time of the estimated auroral electron energy flux (top), 
the height-integrated Joule heating rate (bottom), and the electric potential (bottom, contours) above 50° 
magnetic latitude in the northern hemisphere for a quiet period on 1997 January 10: 0255 UT (left), and an 
active period at 1105 UT (right) on the same day. The electric-potential contour interval is 10 kV (Richmond 
and Lu 2000) 

surges of 100-150 m/s which are predominantly on the nightside. Note further that by this 
time the winds have reached the equator where they will interact with the opposing winds 
from the other pole, producing an apparent “sloshing” of the winds between hemispheres, 
as discussed by these authors. 

Large neutral winds driven by storms can subsequently lead to changes in the plasma 
density, electrodynamics, and composition. Indeed, Lu et al. (2009) showed that disturbance 
neutral winds driven by Joule heating associated with magnetic storms transfer significant 
energy into the mid and low latitudes in the form of gravity waves. These authors showed 
that such gravity waves travel with phase speeds close to the local acoustic speed (~500 m/s 
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Fig. 63 Model neutral wind 
variations driven by auroral 
heating (Fuller-Rowell et al. 
2008). The model calculations 
show the northern component of 
the wind variation in color, 
starting at 3 hours after the 
initiation of an impulse auroral 
heating event. The vectors 
(i arrows ) show the total 
horizontal component, including 
the zonal wind component, 
represented according to the 
scale at the upper right 


Model Wind Variations Driven by Auroral Heating 
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Fig. 64 Model current patterns 
in the northern hemisphere 
showing the disturbance dynamo 
6 hours after the initiation of a 
magnetic storm (after Blanc and 
Richmond 1980) 


Total Horizontal Current 


50mA/m 



or higher) and subsequently push plasma up the magnetic field lines, greatly altering the F- 
region plasma density profile along their track, as observed by these authors in data gathered 
at the Millstone Hill and Arecibo incoherent scatter radar observatories. 

Disturbance Dynamo. Another consequence of the auroral heating and associated winds 
directed equatorward from high latitude regions is the creation of new global circulation 
patterns at lower latitudes, as the equatorward winds develop zonal components at mid- 
latitudes due to Coriolis forces. The resulting circulation creates a system of electric fields 
and currents known as the disturbance dynamo (Blanc and Richmond 1980). During active 
conditions, the Sq current patterns discussed earlier (see Fig. 16) actually split into two cells 
in the upper and lower latitudinal portions of the hemisphere as shown in Fig. 64. Driven by 
auroral heating at high latitudes, a second cell is introduced with an “anti-Sq” rotation in the 
upper middle latitudes. Here, strong poleward electric fields driven by strong poleward elec- 
tric fields that produce a westward E x B drift or eastward currents centered near 45 degree 
latitude. At the equator, the onset of perturbed east-west electric fields has been shown by 
Fejer et al. (1983) and Scherliess and Fejer (1997) to be delayed by appropriate time scales, 
including short (1-12 hrs) and long (22-28 hrs) delays, implying a thermospheric circulation 
source and supporting the disturbance dynamo theory. Other electrodynamic contributions 
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Fig. 65 Images of the earth at 
130.4 nm during the daytime as 
observed by the VIS camera on 
the Polar satellite, showing the 
O/N 2 ratios during a quiet period 
before a storm on April 16, 2002 
(left), and during a magnetic 
storm on April 19, 2002 (right). 
The data are provided by 
J. Sigwarth and J. Kozyra 


Polar Satellite — VIS Earth Camera 


16 Apr 2002 (02/106) 19 Apr 2002 (02/109) 

17:48:24 UT 130.4 nm 19:07:18 UT 130.4 nm 



resulting from magnetic storms include those from penetration electric fields and are dis- 
cussed further on below. 

Composition Changes Related to Magnetic Storms. Other major effects of magnetic 
storms that are quite spectacular are those associated with neutral composition changes that, 
in turn, affect the plasma density. Variations in the column-integrated O/N 2 ratio observed 
by the VIS camera on the Polar satellite are presented in Fig. 65 (J. Sigwarth and J. Kozyra, 
personal communication). The figures show global views of the 130.4 nm airglow on April 
16, 2002 prior to the magnetic storm (left) and on April 19, 2002, during the storm interval 
(right). Such large-scale composition variations are a well-known consequence of magnetic 
storms (e.g., Rishbeth et al. 1987; Prolss 1997; Strickland et al. 2001) and are associated 
with high latitude heating of neutral gases that are then transferred to lower latitudes. The 
heating causes upwelling of molecules in the air (O 2 , N 2 ) and drives changes in the diffu- 
sive equilibrium of the neutral species, enriching heavy species and depleting lighter ones. 
The molecular species become ionized via charge exchange and then undergo rapid loss 
as recombination is much faster for the molecular ions. Ultimately, the lower conductivity 
affects the electrical structure of the ionosphere, as well as its coupling with the magneto- 
sphere. Other examples of changes in the column-integrated O/Nt ratio have been observed 
in the GUVI instrument on the TIMED satellite, as shown in Fig. 66, and compare well with 
simulations using the TIMEGCM model (e.g., Meier et al. 2005). A review of current un- 
derstanding of the variations of O/N 2 ratios during magnetic storms is provided by Crowley 
and Meier (2008). 

Sub-auroral Polarization Streams. Another consequence of magnetic storms is that they 
create increased pressure gradients at the inner edge of the plasma sheet, as well as strong 
conductivity gradients in the ionosphere. This combination results in enhanced, poleward- 
directed electric fields in the pre-midnight sector at the equatorward edge of the auroral oval 
(Wygant et al. 1998). These electric fields map into the ionosphere and can become quite 
large, driving strong westward plasma drifts, particularly in the low conductivity region 
equatorward of the auroral oval in the pre-midnight sector. This phenomenon is known as 
sub-auroral polarization streams or SAPS, a term proposed by Foster and Burke (2002), 
that includes similar phenomena known as polarization jets (Galperin et al. 1973) and sub- 
auroral ion drifts or SAIDS (Spiro et al. 1979) which have been discussed by numerous 
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Fig. 66 Observations of O/No ratios as observed by GUVI on the TIMED satellite for five days from Novem- 
ber 18-22, 2003 (upper panel), and TIEGCM model calculations of this ratio (Meier et al. 2005). The data 
are essentially recorded at the same local time, in this case near noon. Note that the observing time (i.e., UT) 
runs from right to left within each day panel and the panel boundaries are near 0000 UT 


authors (e.g., Smiddy et al. 1977; Anderson et al. 1991, 1993; Karlsson et al. 1998). Situated 
equatorward of the auroral precipitation (hence the term, “sub-auroral”), these streams are 
closely related to where region 2 current are located. In fact, the electric fields are believed 
to increase within the low conductivity regions to maintain Pedersen current closure of the 
field-aligned currents. The appearance of enhanced, latitudinally narrow SAIDS has been 
shown to be related to substorms (e.g., Anderson et al. 1991, 1993; Karlsson et al. 1998). 
The SAPS phenomenon is an important element in the electrodynamics of magnetic storms 
(e.g., Burke et al. 2000). 

Examples of the enhanced electric fields associated with the sub-auroral polarization 
streams measured by the Dynamics Explorer-2 satellite are shown in Fig. 67 together with 
simultaneous measurements of the plasma density (Pfaff et al. 1988). The upper three rows 
show examples of large amplitude (some greater than 300 mV/m) enhancements with narrow 
latitudinal extent, whereas the bottom row shows examples of broader regions of enhanced 
electric fields with somewhat lower amplitudes. Examples of polarization streams that are 
wider and not peaked have also been reported by Foster and Vo (2002). Notice that the en- 
hanced electric fields are often coincident with sharp decreases in the plasma density, though 
this is not always the case. The enhanced electric fields have been observed in conjugate 
hemispheres, as indicated by the pairs of data gathered on the same orbit. These observa- 
tions are not magnetically conjugate as the satellite orbit did not intersect the same magnetic 
field line as it encountered the SAPS event in the northern and southern hemisphere. How- 
ever, the data clearly show that the electric fields are enhanced over a wide range of magnetic 
local time. Simultaneous measurements of this phenomenon at high (9000 km) and low alti- 
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Fig. 67 Examples of large “sub-auroral” electric fields and concurrent plasma density observed by the Dy- 
namics-Explorer-2 satellite. These electric fields are oriented toward the pole producing westward plasma 
drifts often referred to as SAPS (sub-auroral plasma streams). Two pairs of measurements shown here were 
recorded along the same orbit, and reveal similar electric fields in conjugate hemispheres 


tudes (840 km) by probes on the Akebono and DMSP satellites, respectively, were reported 
by Anderson et al. (2001), and support the fact that these fields map between hemispheres 
along the magnetic field. 

SAPS have profound effects on the local ionospheric conductivity and composition, as 
demonstrated in a series of papers by Anderson et al. (1991, 1993). Pointing out that the rate 
of the charge exchange reaction (discussed in Sect. 2): 
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0+ + N 2 -> N0 + + N 

depends on the 4th power of the velocity differential between the ions and the neutrals 
(Schunk and Banks 1976), these authors showed that the N0 + concentration will increase 
substantially in regions where strong SAPS are present. Since the recombination rate of 
N0 + is much faster than 0 + , the plasma number density will decrease further in regions 
of large electric fields. Anderson et al. (1991) also showed that large vertical ion flows that 
result from heating via ion-neutral collisions within SAPS will also contribute to the density 
depletions in the sub- auroral ionosphere. 

6.2 Prompt “Global” Magnetic Storm Effects 

Penetration Electric Fields and Their Effect on Mid- and Low Latitude Plasma Density. The 
broadening of the auroral oval and the transfer of the Joule heating energy to mid-latitudes 
is only part of the story of the upper atmosphere’s response to magnetic storms. For the 
largest storms, the entire ionosphere may be subjected to global scale variations such as 
those driven by prompt penetration electric fields. These have profound effects on the entire 
ionosphere/upper atmosphere at all latitudes. Such penetration fields are associated with 
southward turnings of the IMF and are believed to represent interplanetary electric fields that 
have penetrated the inner magnetosphere to the mid- and low-latitude ionosphere. Examples 
of such fields are shown in Figs. 68a, 68b in which the zonal electric field measured at 
the Jicamarca Radio Observatory in Peru (at the magnetic equator) increases in concert 
with the southward turning of the IMF B z component (Huang et al. 2005). (Note that these 
authors have included a time shift to allow for the solar wind to travel from the satellite 
measurement location in the solar wind upstream of the magnetosphere to the ionosphere, 
which obscures a detailed comparison of the timing of the solar wind events and the exact 
onset of the ionospheric effects.) Although the details of such penetration fields as well 
as the shielding efficiency and time constants of the inner magnetosphere are still being 
debated, these fields appear to have dramatic effects on the distribution of plasma in the low 
and mid-latitudes. 

One of the most important responses of the global ionosphere to large magnetic storms is 
that the low latitude ionospheric plasma can become depleted within a significant latitudinal 
extent centered on the magnetic equator, at least for localized longitude regions. Plasma 
density measurements from consecutive passes of the DMSP satellite in a circular orbit 
at 840 km clearly show how the ionospheric plasma “disappeared” from —20° to +20° 
(magnetic latitude) during a large magnetic storm (at least for the longitudes sampled), as 
shown in Fig. 69 (Greenspan et al. 1991). The simultaneously measured upwards plasma 
drifts were in excess of 200 m/s and are shown by these authors to be consistent with plasma 
being elevated above the satellite altitude. In this manner the phenomenon is similar to the 
fountain effect (discussed in Sect. 4), driven by sustained, large amplitude east-west zonal 
electric fields. Discussed in terms of a so-called “super fountain” effect, Heelis (2008) argues 
that, taking into account neutral winds and other effects, most of the plasma transport under 
such conditions is poleward (i.e., horizontal). When such density “holes” are created during 
the daytime due to the large plasma uplift, note that photo-ionization continues at the lower 
altitudes and the creation of a new F-region proceeds. 

An example of large increases in the plasma density in the upper F-region during the 
daytime occurring within a significant magnetic storm is provided by measurements of total 
electron content (TEC) above 400 km, as measured by the CHAMP satellite and shown in 
the lower panel of Fig. 70 (Mannucci et al. 2005). Here, data from three consecutive satel- 
lite orbits (spaced ~90 minutes apart) show how the TEC changes during the storm main 
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Solar Wind Observations -- ACE Satellite 
Ionospheric Electric Fields -- Jicamarca Obs., Peru 


R.F. Pfaff 


January 13, 1999, event day 
January 12, 1999, quiet day 


November9-10, 2004, event day 
November 11-12, 2004, quiet day 



UT (Hours) 
(a) 


UT (Hours) 
(b) 


Fig. 68 Solar wind IMF B z (top) measured by the ACE satellite with associated interplanetary electric 
field, Ey ( middle ) and the equatorial zonal ionospheric electric field (IEF) E zona \ measured by the Jicamarca 
radar (bottom). Data from two storm periods are shown, for January 13, 1999 (left) and November 9-10, 
2004 (right). Also shown are data recorded during quiet conditions on neighboring days for comparison. The 
periods shaded in red and blue highlight the interval of positive IEF, corresponding to southward IMF, and 
enhanced eastward ionospheric electric fields, respectively. Notice the increased scales in the panels on the 
right, both for the solar wind measurements as well as for the electric fields measured at Jicamarca (adapted 
from Huang et al. 2005) 


phase, indicated by the decrease in Dst, which is plotted in the upper panel. Notice how the 
maxima of the equatorial anomaly that are normally observed nearer to ±15° (see blue trace 
corresponding to 19:00 UT) were “pushed” to values of ±30° (see red and black traces) yet 
still retained their approximately symmetric, winged shape. These data were gathered near 
1300 hrs SLT at low and middle latitude (A. Mannucci, personal communication, 2012) 
and show prompt dayside total electron content increases of ~900 % above 400 km, sug- 
gesting that the normal, daytime ionosphere at lower altitudes had risen above the CHAMP 
satellite altitude. These authors also invoke prompt penetrating electric fields which create 
plasma uplift with associated diffusion along the magnetic field lines to higher latitudes to 
explain both the large TEC values and the observed, symmetric latitudinal spreading. A dif- 
ferent interpretation for the storm-time density enhancements on the dayside is that due to 
high latitude convection electric fields expanded equatorward, as discussed by Heelis et al. 
(2009a, 2009b) and Heelis and Mohapatra (2009). 

Examples of Low and Mid-Latitude Irregularities Associated with Magnetic Storms. To 
demonstrate further how the ionospheric storm response appears to be initiated at the equator 
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Fig. 69 Log of total ion density 
in cm~^ ( upper panels, light 
orange ) measured by the DMSP 
spacecraft plotted as a function of 
magnetic latitude for live 
successive passes, each near 
2130 LT, during the magnetic 
storm of March 13-14, 1989. The 
lower panels show measurements 
of the vertical ion drifts in which 
the positive values are shaded 
blue. The labels to the right of 
each plot give the universal time 
and geographic latitude and 
longitude at which the spacecraft 
crossed the magnetic equator 
(adapted from Greenspan et al. 
1991) 


DMSP F9 NIGHT 



March 13 
21:20 UT 
1 3° GLAT 
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March 13 
23:03 UT 
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March 14 
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March 14 
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Magnetic "Superstorm" of Oct. 30-31, 2003 



Dayside TEC Above CHAMP Satellite at 400 km 



-60 -40 -20 0 20 40 60 

Magnetic Latitude (Dipole) After Mannucci et al, 2005. 


Fig. 70 The integrated electron content (corresponding to total electron content or TEC) as measured by the 
GPS instrument on the CHAMP satellite for altitudes above the CHAMP altitude of 400 km ( lower panel). 
Data are shown for three orbits plotted as a function of magnetic latitude which correspond to times just prior 
to (blue) and after (red and black) the onset of a magnetic storm on October 30, 2003, as shown in the plot of 
Dst in the upper panel. The locations of the CHAMP orbits are shown in the upper right corner of the lower 
panel. The local times of these orbits range from 1230-1330 LT for latitudes within ±60 degrees. Points 
missing near the anomaly trough are due to the elevation angle cut-off (adapted from Mannucci et al. 2005) 


and then move polewards simultaneously in both hemispheres, at least during the evening, 
we show DEMETER satellite data in Fig. 7 1 . Consecutive orbits of plasma density and elec- 
tric field irregularities are presented for two different magnetic storms, as shown in the Dst 
data in the upper panels (Pfaff et al. 2008). These data were gathered along circular orbits at 
7 10 km in the evening (~22.5 LT) sector. At the storm commencement the first indication of 
both enhanced plasma density and associated irregularities at 710 km occurred at the equa- 
tor, and then proceeded towards the poles in both hemispheres as the storm progressed. The 
irregularities associated with the inter-hemispheric, unstable flux tubes continue to rise in 
altitude and latitude, symmetrically in both hemispheres and are encountered by the consec- 
utive DEMETER orbits at successively higher latitudes. These data reinforce the hypothesis 
that the storm response at low and mid-latitudes is initiated at the equator. At the same time, 
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Fig. 71 DEMETER satellite data (near 22:30 LT) for the magnetic storm periods of August 24, 2005 (a) and 
May 15, 2005 (b). The top panels show the Dst estimates based on ground-based magnetic field measure- 
ments. The middle and lower panels show the DEMETER integrated electric field irregularity power from 
10-100 Hz and the plasma number density along each orbit. The satellite data have been converted to color 
scales and mapped to bins of magnetic latitude between —65 and +65 degrees (Pfaff et al. 2008) 


the polar cap high latitude ionosphere widens and expands equatorward, as evidenced in 
the lower panels in Fig. 7 1 where the plasma density decrease or trough near 60 degrees 
magnetic latitude prior to the storms expands equatorward towards ±40 degrees magnetic 
latitude in both the northern and southern hemispheres. 

Ground-based observatories at mid latitudes have also observed ionospheric irregular- 
ities associated with magnetic storms that provide important clues concerning whether 
unstable flux tubes might be set up between hemispheres, as suggested by the DEME- 
TER data above. Two examples of radar measurements gathered in Puerto Rico during 
magnetic storms are shown in Fig. 72. The upper right-hand panel shows data gathered 
with a VHF coherent scatter radar looking northward from Puerto Rico (figure cour- 
tesy E. Kudeki; see Swartz et al. 2000) and shows irregularities emerging from what 
is believed to be descending F-region ledge, similar to those associated with low lat- 
itude spread-F. The lower panel shows Arecibo incoherent scatter measurements dur- 
ing a similar magnetic storm (Kelley et al. 2002). Although the fixed radar beam can- 
not distinguish temporal and spatial variations that may be at work during this event, 
the data suggest the presence of depleted flux tubes that emerge from the F-region 
ledge and which may be expected to give rise to electrostatic irregularities. Other mea- 
surements (e.g., TEC and airglow depletions) during these and similar magnetically- 
disturbed events at mid latitudes have been reported by other researchers (e.g., Makela 
et al. 2000; Kelley et al. 2000; and Bust et al. 2000) and demonstrate that the mid- 
latitude ionosphere becomes highly structured at night during magnetically disturbed con- 
ditions. 
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VHF Backscatter Radar Power 
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(b) Geophysical Indices -- June 24-28, 1998 



Electron Density 

Arecibo Observatory, Puerto Rico 



0000 0100 0200 0300 0400 



June 26, 1998 -- Local Time (AST) 


KeBey etal. [2002] 


Fig. 72 Measurements by ground-based instruments in Puerto Rico of plasma density structures associated 
with two modest magnetic storms, whose Dst and Kp signatures are shown on the left hand side. The panel 
on the upper right shows 50 MHz backscatter radar measurements that resemble spread-F scatter that were 
coincident with the period of largest negative Dst (figure courtesy E. Kudeki; See also Swartz et al. 2000). The 
panel on the lower right shows plasma density measurements from the incoherent scatter radar at Arecibo 
that were also coincident with the period of the largest negative Dst (Kelley et al. 2002). The large plasma 
depletion “chimneys” that emerge from the lower edge of the F-region between 0100 and 0200 LT resemble 
spread-F depletions more commonly observed at lower latitudes 


The panels showing the Dst data on the left hand side of Fig. 72 indicate that in each 
case, the radar irregularities and depletions were observed at the peak times of their respec- 
tive storms. The fact that the radars in each case were situated on the nightside during each 
of these two storm events, in which the main phase lasted less than 12 hours, undoubtedly 
helped enable the observations of these particular storm signatures. These mid-latitude ob- 
servations of F-region irregularities and depletions resemble spread-F processes known to 
occur at low latitudes at night. 

Storm-Enhanced Density Structures. A final example of significant effects on the near- 
earth plasma from magnetic storms is planetary-scale plasma density structures that ap- 
pear to break off from the enhanced mid-latitude plasma that accumulates during a 
large amplitude magnetic storm and then is transported poleward, as shown in Fig. 73 
(figure courtesy of A. Coster). These plumes are referred to as storm enhanced den- 
sity (SED) erosion plumes (see Foster et al. 2005, and references therein) and are evi- 


•£) Springer 





The Near-Earth Plasma Environment 


101 


Fig. 73 Storm-enhanced plasma 
density (SED) signatures in total 
electron content (TEC) observed 
on November 20, 2003. These are 
believed to be connected to 



20 November 2003 
19:30:00 to 19:45:00 U.T. 


TEC 


50 


plasmasphere erosion and driven 
by sub-auroral electric fields 


from the inner magnetosphere. 


Strong plasma density gradients 
are observed by a network of 


ground-based GPS observatories. 
Figure courtesy of A. Coster 


dent in widely spaced GPS measurements of total electron content. Foster et al. (2005) 
reported that the plume shown in the figure is carried through the cusp during the 
main phase of the magnetic storm of November 20, 2003. In such a scenario, the en- 
hanced plasma density is driven upwards and poleward by the prompt penetration elec- 
tric fields, then towards the dayside by storm-related SAPS, and finally by dayside, anti- 
sunward convection to form the plume across the polar cap at higher latitudes. Such 
plumes which arrive at the high latitudes may enhance ion outflow into the magneto- 
sphere. 

7 Summary 

In this paper, we have reviewed the origins of the near-earth plasma created by photo- 
ionization due to solar EUV radiation, primarily at mid and low latitudes, and by impact 
ionization from precipitating energetic particles at higher latitudes. Not only is the near- 
earth plasma created directly from the earth’s extended atmosphere, but the ionization that 
is formed remains coupled to the neutral gases and their dynamics. Clearly, to fully under- 
stand the near-earth space environment, we must consider both the neutral and ionized gases 
as an ensemble, forming one envelope that surrounds our planet. 

We have examined many of the numerous processes and dynamics inherent in the cou- 
pled ion-neutral gas at low and mid latitudes. Throughout this region, the ionized gas com- 
ponent is largely controlled by the neutral dynamics. However, as we have shown, small 
changes between the ion and electron distributions, primarily set up as a result of collisions 
with the neutral particles, are sufficient to generate global systems of currents and elec- 
tric fields. These currents map along the magnetic field lines to other regions of Geospace 
where their effects may be profound, for example near the magnetic equator just after sun- 
set. Furthermore, the fact that the ionosphere/upper atmosphere is strongly influenced from 
forcing from below — i.e., from tides, gravity waves, and thunderstorm activity originating 
in the troposphere — underscores the controlling role of the atmosphere on the behavior of 
the ionized gases, particularly at mid and low latitudes. 
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We have also examined the dynamics and processes inherent in the coupled ion-neutral 
gas at high latitudes. We have shown how strong electric fields and currents originating in 
the magnetosphere often dominate the dynamics, “stirring the pot” of the polar regions with 
complex convection and field-aligned current patterns that together yield a self-consistent 
electrodynamic picture. As the magnetospheric electric fields drive neutral winds via ion- 
neutral collisions, the neutral atmosphere is shown to play a key role in the conversion of 
magnetospheric electromagnetic energy to mechanical energy in the high latitude upper at- 
mosphere. Furthermore, the neutral atmosphere also plays a critical role in determining how 
the magnetospheric energy is dissipated, primarily through Joule heating which essentially 
results from friction between the differential motion of the neutral and ionized gases. In- 
deed, the response of the neutral gases includes large horizontal and vertical flows that are 
not confined by the magnetic field and which can subsequently convey the magnetospheric 
energy to locations far beyond where it was initially deposited. The neutral atmosphere also 
affects the closure, in the lower ionosphere, of magnetospheric field-aligned currents via 
Pedersen currents. Finally, the local conductivity in the ionosphere, which is also affected 
by the neutral atmosphere, has been shown to have a direct effect on the creation of elec- 
trostatic structures at higher altitudes that accelerate the energetic particles that ultimately 
precipitate to form the discrete aurora. 

We have also examined many of the fundamental physical processes inherent in the near- 
earth plasma that result from magnetic storms, such as the effects of increased energy depo- 
sition at high latitudes, and have shown how such energy is then distributed away from the 
polar regions. We also showed evidence that the global ionosphere is subjected to large scale 
variations when the earth is impacted by magnetic storms of appreciable magnitude. Such 
effects include substantial changes in the ionospheric density, which may be attributable to 
prompt penetration electric fields which have been observed at the magnetic equator near 
the onset of magnetic storms. These observations demonstrate the global and rapid response 
of the nearearth plasma to solar wind forcing associated with significant magnetic storms. 

Progress in our understanding of the ionosphere/thermosphere system requires adopting 
a systems approach with respect to interpreting the large scale dynamics associated with the 
energy and momentum exchanges within this gaseous envelope as it adjusts internally and 
responds to input from both the troposphere below and the magnetosphere above. Under- 
standing such a self-consistent electrodynamic and hydrodynamic response to these drivers 
is a critical, long-standing objective which is at the forefront of future research in this field. 

One area that is highly relevant to future progress in our understanding is that of the 
tremendous promise within the modeling community of incorporating both physics and 
observations in computational models of the upper atmosphere. Examples of such mod- 
els that include first principles of ionosphere/thermosphere coupling are the Thermosphere - 
Ionosphere-Mesosphere-Electrodynamics General Circulation Model (TIMEGCM) (Roble 
et al. 1988) and the Coupled Thermosphere-Ionosphere-Plasmasphere-Electrodynamics 
(CTIPe) model (Fuller-Rowell et al. 1996). Furthermore, the Whole Atmosphere Model 
(WAM) (Fuller-Rowell et al. 2010) is an extension of the US National Weather Service 
Global Forecast System and includes the impact of tides, planetary waves, and gravity 
waves. Finally, the Assimilative Mapping of Ionospheric Electrodynamics (AMIE) model 
(Richmond et al. 1992) includes the most accurate high latitude specification, which is im- 
portant for understanding the high latitude forcing. 

While the models mentioned above rely on both physical understanding and observa- 
tions, in some areas the models are considerably “data starved”. For example, very few 
simultaneous measurements of the dynamics of neutral and ionized gases at the lower alti- 
tudes exist. Although targeted research rockets, satellites, and ground-based measurements 
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are essential for advancing our understanding, the systems approach requires simultaneous 
measurements using a number of observational techniques, including space-based in situ 
and imaging experiments as well as arrays of ground-based radars and other instruments. 
The space-based direct measurements show particular promise when organized as “constel- 
lation” missions, particularly in configurations where data may be gathered simultaneously 
at multiple local times and longitudes. Such platforms promise to reveal how the neutral 
and ionized dynamics are organized, as well as show how energy and momentum is dis- 
tributed, and re-distributed, globally. For example, such measurements would reveal the 
two-dimensional coupling of neutral and plasma gases and their dynamics across the entire 
high latitude ionosphere, show the global consequences of propagating tides, and demon- 
strate how magnetic storms promptly effect the entire ionosphere/thermosphere system as a 
function of latitude, local time, and longitude. 

Finally, we emphasize that to be most effective, such measurement platforms need to 
journey where the most important physical processes are actually taking place. In other 
words, to examine how the ionized and neutral gases exchange energy and momentum, it is 
important that the instruments sample altitudes where these interactions are most profound. 
For example, simultaneous measurements of ions and neutrals below 400 km are essential 
to understand the coupled behavior of these gases. Focused missions that furthermore ex- 
plore well below 300 km are expected to reveal important fundamental physical processes 
involving the ion-neutral gases where the coupling is maximized and where our knowledge 
of the resulting gas dynamics is rudimentary, at best. 

Ultimately, the behavior of the Earth’s ionized and neutral gases depends on non-linear 
coupling between complex, highly varying external sources. These external sources include 
the variable EUV radiation from the sun as well as intricate energy and momentum ex- 
changes with the troposphere/stratosphere/mesosphere below and the magnetosphere and 
solar wind above. It is precisely because these coupled regions are exceedingly rich in phys- 
ical processes that the near-earth plasma is one of the most exciting, and important, natu- 
ral environments in which to focus scientific research and resources. The new knowledge 
promised from such future research is essential, not only for understanding the Earth, but 
also for understanding the upper atmospheres of other planets, both in our solar system and 
beyond. 
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